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Under  contract  N0(X)14-89-C-2265,  SFA  provided  support  to  the  Naval  Research  Laboratory  Beam  Physics 
Branch  by  researching  and  developing  Uie  Modified  Betatron.  This  branch  is  engaged  in  research  and 
development  of  compact,  high-power  accelerators  and  their  aj^lications.  The  major  research  emphasis  is  on 
high-current,  cyclic  accelerators,  such  as  the  modified  betatron  located  at  NRL  in  Washington,  DC.  In  a 
betatron,  electrons  are  accelerated  to  high  energies  by  an  electric  field.  Particular  applications  of  interest 
include  intense  coherent  radiation  sources  and  beam  propagation. 

SFA  has  supported  the  Beam  Physics  Branch  for  12  years.  Under  multiple  tasks,  SFA  technical  personnel 
have  engineered,  designed,  fabricated,  installed,  and  tested  special-purpose  equipment,  components,  and 
instruments  for  conducting  various  research  projects  to  improve  the  capabilities  of  compact  accelerators. 
This  engineering  and  development  work  is  performed  in  close  collaboration  with  the  experimental  effort 
to  ensure  compatibility,  optimal  function,  and  effective  use  of  government  equipment  and  facilities. 

The  goal  of  the  modified  betatron  accelerator  project  was  to  accelerate  a  1  kA  electron  beam  to  an  energy  of  50 
Me V.  The  research  program  centered  on  the  critical  physics  issues  associated  with  this  goal,  such  as 
beam  stability,  self  field  effects,  image  forces,  and  injection  and  extraction  processes.  Achievement  of  this 
goal  revolutionized  the  DoD’s  directed-energy  weapons  technologies,  establishing  the  Beam  Physics 
Bnnrn’s  modified  betatron  actx;Ierator  as  the  leehi illogical  leader  in  the  field.  SFA  has  played  an 
integral  supporting  role  on  this  project  from  the  novel  construction  of  the  vacuum  chamber  to  expert 
operation  and  day-to-day  maintenance  and  modification  of  the  betatron. 

The  modified  betatron  accelerator  is  a  large,  complex  facility.  Several  high-energy  storage-capacitor  banks 
and  pulsed-power  systems,  with  the  associated  distribution  and  coil  systems,  are  required  to  create  the  various 
magnetic  fields  for  acceleration  and  beam  confinement  and  to  supply  the  beam.  In  addition  to  the  electrical 
components,  the  structural  framework  to  support  the  betatron  chamber  and  associated  magnetic  field  coils 
must  allow  for  precision  positioning  of  the  various  coils  and  yet  be  able  to  withstand  considerable  forces  and 
torques  during  a  shot.  Due  to  this  complexity,  extensive  mechanical  and  electrical  engineering  was  required 
for  the  modification  and/or  fabrication  of  new  parts  for  the  betatron.  SFA’s  engineering  support  consisted  of: 

•  Design  Engineering-ln  cooperation  with  government  scientists,  SFA  provided  support  to  conceptualize 
research  equipment  that  satisfied  the  structural,  mechanical,  electrical,  vacuum,  magnetic,  electrostatic, 
cost,  and  schedule  requirements  of  the  project.  SFA  prepared  layout  drawings,  performed  analyses  for 
such  properties  as  strength  and  deflection,  selected  materials,  and  prepared  detail  drawings. 

•  Developmeni-W^en  required,  SFA  undertook  the  development  of  specific  technologies  that 
were  necessary  to  successfully  build  required  research  equipment.  SFA  identified  materials  and 
located  sources,  located  fabricators,  and  built  prototypes.  SFA  obtained  competitive  bids  and 
placed  purchase  requisitions,  monitored  the  progress  of  fabricators,  and  took  delivery  of  equipment 
as  required.  SFA  also  transported  equipment  to  NRL;  inspected  the  equipment;  and  assisted 
in  installation,  preliminary  adjustments,  and  operation.  This  area  also  included  fabricator  development 
In  the  past,  research  equipment  fabrication  requirements  arose  for  which  there  was  no  known  industrial 
base  or  the  existing  base  was  too  expensive.  SFA,  therefore,  worked  with  various  contractors  or 
subcontractors  as  required  to  assist  in  the  development  of  the  technology  base.  This  x'as  accomplished 
by  locating  contractors/subcontractors  that  possessed  related  skills,  then  educating  them  to  develop  the 
required  base. 

•  Testing-SFA  undertook  extensive  testing  of  individual  components  and  materials  to  define  and  verify 
important  characteristics  relevant  to  the  proposed  use.  SFA  also  maintained  databases  of  test  results, 
including  both  in-house  testing  and  testing  carried  out  at  other  sites.  The  development  of  new  or 
modification  of  existing  systems  entails  the  design,  assembly,  testing,  and  installation  of  new 
components  or  modifications  to  existing  components  to  improve  the  performance  of  the  accelerator  and 
to  investigate  key  physics  issues  (such  as  extraction).  Branch  scientists  were  responsible  for  the 
conceptual  design  of  equipment.  Through  interaction  with  SFA's  engineers,  the  conceptual  design  was 
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refined  and  detailed  designs  produced.  During  this  process,  manufacturability  of  the  equipment  was  also 
addressed.  In  particular,  a  straightforward  design  process  can  result  in  difficult-to-manufacture 
components.  To  ensure  that  designs  were  producible,  SFA  investigated  several  possible  manufacturing 
processes  at  the  outset.  Using  this  approach,  SFA  often  arrived  at  novel  manufacturing  solutions.  An 
example  of  this  was  SFA’s  participation  in  the  MBA  vacuum  chamber  development.  Several  materials 
were  considered  for  the  chamber,  and  SFA  engineers  determined  that  a  fibetgl ass/epoxy  chamber  would 
meet  the  project  needs.  To  manufacture  the  components,  a  local  fiberglass  boat  manufacturer  was  trained 
by  SFA  to  produce  the  toroidal  vacuum  chamber  sections.  This  effort  resulted  not  only  in  a  superior 
product  but  also  in  a  reduction  in  cost. 

In  support  of  this  contract,  SFA  has  performed  numerous  tasks  related  to  vacuum  systems.  These  tasks  have 
included  vacuum  chamber  design,  leak  checking,  and  measurement  of  outgassing  rates  and  pumpdown 
rates.  Recently.  SFA  personnel  designed  a  cylindrical  vacuum  chamber  with  a  pumped  volume  of  0.5  m3 
out  of  304  stainless  steel  alloy.  Thirty  ports  were  designed  to  allow  for  electron  gun  introduction, 
optical  observation,  pumping  connections,  mechanical  manipulation,  and  various  diagnostics.  These  ports 
ranged  in  diameter  from  1  inch  to  8  inches.  Single-side  vacuum  welds  were  specified  to  eliminate  the 
possibility  of  virtual  leaks-trapped  pockets  of  gas  on  the  interior  of  a  vacuum  system.  Stress  and  stability 
analyses  were  carried  out  to  ensure  structural  integrity.  Procurement  drawings  were  prepared  and  competitive 
bids  were  obtained.  A  contract  was  placed  with  a  reasonable  price  and  delivery  schedule.  During 
fabrication  of  the  vacuum  chamber,  design  work  was  completed  on  the  suppon  structures,  baking  and 
cooling  systems,  and  pump  selection.  After  delivery  of  the  vacuum  chamber,  the  system  was  assembled, 
groomed,  and  placed  into  operation.  A  single  Sargent- Welch  Turbo  Molecular  Pump  Station  having  a 
capacity  of  400 1/s  was  able  to  pump  the  system  down  to  2  x  1 0'^  torr  in  approximately  4  hours.  With  baking 
and  further  grooming,  a  vacuum  of  7  x  10  ®  torr  was  achieved. 
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Bean  Extraction  Scheme  from  the  Modified  Betatron  Accelerator 

C.  A.  Kapetanakos,  S.  J.  Marsh,^^  and  D.  Diaietis^'’^ 

Plasma  Physics  Division,  Naval  Research  Laboratory,  fyashinpon.  D.C.  20375 
(Received  10  March  1988) 

A  technique  is  proposed  for  extraction  of  the  electron  ring  from  the  modiSed  betatron  accelerator. 
Basically,  this  technique  consists  of  exciting  the  resonance  that  naturally  exists  for  some  speciSc  values 
of  the  ratio  of  the  vertical  to  toroidal  magnetic  field. 


PACS  numbers:  S2.75.Di.  29.20.Fj,  41.80.E« 

The  modified  betatron  accelerator'-^  is  one  among  the 
several  compact,  high-current  accelerator  concepts  cur¬ 
rently  under  development  in  various  laboratories.  In  this 
device  a  strong  toroidal  magnetic  field  Bg  has  been  added 
to  the  conventional  betatron^  magnetic  field  configura¬ 
tion.  Although  Bg  substantially  improves  the  stability  of 
the  conventional  betatron,  the  beam  injection  and  cap¬ 
ture  and  the  electron  ring  extraction  after  the  completion 
of  acceleration  are  substantially  more  involved  as  a  re¬ 
sult  of  the  toroidal  field. 

In  this  Letter,  we  report  on  an  extraction  scheme  that 
is  easily  realizable  and  has  the  potential  to  lead  to  very 
high  extraction  efficiency.  Briefly,  the  proposed  extrac¬ 
tion  scheme  is  based  on  the  transformation  of  the  circu¬ 
lating  electron  ring  into  a  stationary  helix,  in  the  toroidal 
direction,  by  excitation  of  the  resonance  that  naturally 
exists  for  some  specific  values  of  the  ratio  of  the  vertical 
to  toroidal  magnetic  field.  Transformation  of  the  ring 
into  a  helix  is  achieved  with  a  localized  vertical  magnetic 
field  disturbance  that  is  generated  by  an  agitator  coil. 
As  the  minor  radius  of  the  helix  increases  with  each  pas¬ 
sage  through  the  gap  of  the  agitator  coil,  the  electrons 
eventually  reach  the  extractor,  which  has  the  property 
that  all  the  magnetic  field  components  transverse  to  its 
axis  are  equal  to  zero.  Thus,  the  electron  ring  upwinds 
into  a  straight  beam. 

Although  the  proposed  scheme  is  related  to  the  stan¬ 
dard  techniques^-’  used  for  extraction  of  the  beam  from 
existing  circular  accelerators,  several  of  its  key  features 
are  different  because  of  the  presence  of  the  toroidal  mag¬ 
netic  field. 

Extraction  scheme. — After  the  completion  of  accel¬ 
eration,  i.c.,  when  the  desired  electron  beam  energy  has 
been  achieved,  the  electron  ring  centroid  is  displaced  ra¬ 
dially  by  the  intentional  mismatch  of  the  magnetic  flux 
and  the  betatron  magnetic  field.  In  the  results  that  will 
be  shown  in  the  next  section,  this  mismatch  has  been 
achieved  by  superimposition  of  a  low-amplitude  vertical 
magnetic  field  that  varies  exponentially  with  time  on  the 
betatron  field.  It  has  been  shown  theoretically  and 
verified  by  extensive  numerical  results  that  during  the  ra¬ 
dial  displacement  of  the  ring  centroid  the  amplitude  of 
the  slow  mode'  remains  very  small,  i.e.,  a  few  millime¬ 
ters,  provided  that  the  mismatching  field  varies  slowly 


with  respect  to  the  ring  bounce  (poloidal)  period.  Fur¬ 
thermore,  computer  simulations  with  the  U.S.  Naval 
Research  Laboratory  mobe  particle-in-cell  computer 
code  have  shown  that  during  the  radial  displacement, 
that  lasts  several  microseconds,  the  minor  cross  section 
of  the  ring  preserves  its  integrity  and  the  ring  emittance 
remains  constant. 

As  the  major  radius  of  the  ring  centroid  increases 
slowly  with  time,  the  gyrating  electrons  reach  the  local¬ 
ized  magnetic  disturbance  generated  by  the  agitator  coil. 
At  this  radial  position  the  ratio  of  the  vertical  magnetic 
field  Bt  to  the  torodial  magnetic  field  Bg  has  been  select¬ 
ed  to  satisfy  the  condition 

B:/Bg~2in2l^-\),  (1) 

where  /  “1, 2, 3 _ 

Equation  (1)  implies  that  the  frequency  of  the  fast 
mode'  is  /  times  the  frequency  of  gyration  around  the 
major  axis.  When  Eq.  (1)  is  reduced  to 

flg^lClz,  where  (ig^eBgJm  and  Ciz  ^eBjm. 

The  purpose  of  the  magnetic  disturbance  is  to  excite 
the  resonance.*-^  As  an  electron  enters  the  lower  mag¬ 
netic  field  region  of  the  disturbance,  its  velocity  vector, 
which  initially  is  directed  in  the  toroidal  direction,  ro¬ 
tates  slightly  in  the  radial  direction,  i.e.,  the  electron  ob¬ 
tains  a  radial  velocity  component.  It  can  be  shown  from 
the  equations  of  motion  that  this  radial  velocity  is  given 
by 

-2(Ani‘'/y)raA9,  (2) 

where  Anf  is  the  cyclotron  frequency  that  corresponds 
to  the  field  of  the  disturbance  generated  by  the  agitator 
coil,  y  is  the  relativistic  factor,  Ca  is  the  radial  distance  of 
the  agitator  coil,  and  Afl  is  the  toroidal  half  width  of  the 
magnetic  disturbance. 

As  a  result  of  the  acquired  radial  velocity,  the  elec¬ 
trons  start  to  gyrate  in  the  torodial  magnetic  field  with  a 
radius 

p-2(A'//)(An//n.)rflA0,  (3) 

where  N  is  the  number  of  passes  through  the  distur¬ 
bance.  If  condition  (1)  is  not  satisfied,  p  grows  as 
instead  of  proportionally  to  A. 
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Since  y  is  very  large,  self-fields  can  be  ignored.  How¬ 
ever,  because  of  the  gradient  of  Bt  the  slow  mode' 
(bounce  motion)  is  still  excited  and  the  orbits  of  elec¬ 
trons  in  the  transverse  (/■.z)  plane  process  very  slowly. 
Therefore,  for  times  short  in  comparison  with  the  bounce 
period,  i.e.,  for  a  few  revolutions  around  the  major  axis, 
ail  the  electrons  of  the  ring  perform  coherent  motion  and 
a  stationary  helix,  in  the  toroidal  direction,  is  formed.  A 
top  view  of  the  helix  is  shown  in  Fig.  1,  for  /■■3. 

Ideally,  the  radial  gradient  of  the  magnetic  distur¬ 
bance  should  be  extremely  high,  because  otherwise  the 
fast  mode '  is  excited  before  the  ring  reaches  the  distur¬ 
bance.  In  the  computer  runs  of  the  next  section,  a  dis¬ 
turbance  with  a  satisfactorily  sharp  radial  gradient  is  ob¬ 
tained  by  the  single-turn  agitator  coil  shown  in  Fig.  2. 
The  radial  gradient  of  the  disturbance  is  further  im¬ 
proved  with  two  single-turn  loops  that  are  located  at  the 
edges  of  the  gap.  In  the  computer  runs,  the  magnetic 
field  of  the  disturbance  has  been  obtained  from  exact 
analytical  expressions  that  are  too  long  to  be  given  here. 

With  successive  passes  through  the  disturbance  of  the 
agitator  the  radial  excursion  of  the  orbit  increases  until 
the  gyrating  electrons  reach  the  extractor,  which  is  locat¬ 
ed  at  and  at  a  slightly  greater  radial  distance  than 
the  agitator  coil.  The  results  of  the  next  section  were  ob¬ 
tained  with  a  simple  extractor  consisting  of  two  parallel 
plates  with  current  flowing  in  opposite  directions.  These 
two  plates  have  infinite  extent  in  the  z  and  semi-infinite 
extent  in  the  y  directions.  The  linear  current  density  of 
the  plates  is  adjusted  to  make  the  total  Bt  between  the 
plates  at  0“O  equal  to  zero.  The  side  of  the  extractor  at 
0”O  is  completely  enclosed.  As  a  result  the  fringing 
fields  are  absent.  The  electrons  enter  the  extractor 
through  a  thin  conducting  foil.  At  the  entrance  of  the 
extractor  the  vertical  displacement  of  the  electrons  and 
their  radial  velocity  arc  almost  zero.  However,  they  have 
a  small  vertical  velocity. 


FIG.  I.  Schemaiic  of  the  proposed  extraction  scheme. 


In  practice,  this  extractor  can  be  realized  by  our  bend¬ 
ing  the  two  plates  to  form  a  torus.  In  order  for  the  field 
to  be  uniform  over  a  finite  vertical  distance,  the  cross 
section  of  each  plate,  after  bending,  should  be  D  shaped. 
In  the  results  of  the  next  section,  the  orbit  of  the  extract¬ 
ed  beam  is  terminated  after  it  propagates  tens  of  centi¬ 
meters  inside  the  extractor.  The  reason  is  that  is  in¬ 
dependent  of  y  while  the  betatron  field  decreases  with  y. 
Thus,  cancellation  of  the  fields  is  not  achieved  over  the 
entire  length  of  the  extractor.  In  praaice  exact  cancella¬ 
tion  of  the  two  fields  can  be  obtained  by  an  increase  in 
the  separation  of  the  two  plates  as  y  increases. 

In  the  previous  discussions,  we  have  assumed  (hat  the 
magnetic  disturbance  generated  by  the  agitator  coil  is 
static.  An  alternative  mode  of  operation  is  to  expand  the 
ring  until  it  reaches  the  gap  of  the  agitator  coil  and  then 
to  pulse  the  coil  rapidly.  Since  the  inductance  of  the  agi¬ 
tator  is  typically  only  a  few  nanohenries,  short  rise  times, 
of  the  order  of  a  few  nanoseconds,  can  be  achieved  with 
modest  voltages.  In  the  pulsed  mode  of  operation  the 
fraction  of  the  ring  that  will  be  lost  is  approximately 
equal  to  the  ratio  of  coil  rise  time  to  period  of  gyration 
around  the  major  axis. 

Finally,  it  should  be  noticed  that  an  ion  channel^ 
formed  by  a  laser  beam  along  the  axis  of  the  extractor 
may  improve  the  extraction  process  and  eliminate  the 
need  for  an  additional  coil  to  cancel  the  component  of  B0 
that  is  transverse  to  the  axis  of  the  extractor  or  the  need 
to  cancel  completely  the  B:  inside  the  extractor. 

Results.— Vi/e  have  studied  the  proposed  extraction 
scheme  in  both  the  static  and  pulsed  modes  for  a  range 
of  parameters  that  arc  compatible  with  the  U.S.  Naval 
Research  Laboratory  modified  betatron  accelerator.  In 
this  Letter,  we  present  results  from  three  runs:  one  in 
the  pulsed  mode  for  y"'40  and  two  in  the  static  mode  for 
y"40  and  400.  The  various  parameters  of  these  three 
runs  are  listed  in  Table  I.  Since  y'»  1,  self-fields  and  im- 


Agltator  Col) 

FIG.  2.  Agitator  coil  that  generates  the  localized  distur¬ 
bance.  It  IS  powered  by  a  coaxial  transmission  line 
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TABLE  I.  Lilt  of  various  parameters  for  the  runs  shown  in  Figs.  3  to  5. 


Run  No. 

Parameter 

267 

266 

268 

Agitator's  mode 

Pulsed 

Static 

Sutic 

Relativistic  factor  y 

40 

40 

400 

Major  radius  ro  (cm) 

100 

100 

100 

Vertical  field  at  ro  (G) 

649.9 

649.9 

6S0I 

Toroidal  field  at  ro  (G) 

-1921 

-1971 

-19700 

Field  index  n 

0.5 

0.5 

0.5 

Resonance  integer  / 

3 

3 

3 

Amplitude  of  mismatching  field  (G) 

— 

60 

600 

Time  constant  of  mismatching  field  (psec) 

.  .  . 

10 

10 

Agitator 

Toroidal  position 

1.3)r 

1.26* 

1.26* 

Toroidal  width  2AS  (rad) 

0.05 

0.066 

0.066 

Inner  radius  (cm) 

120 

120 

120 

Outer  radius  (cm) 

122 

124 

124 

Opening  (cm) 

1.0 

2 

2 

Linear  current  density  (kA/cm) 

0.25 

0.375 

3.750 

Field  ABf  (G) 

-300 

-450 

-4500 

Extractor 

Opening  toroidal  position 

0 

0 

0 

Minimum  inner  radius  (cm) 

121.5 

120.5 

120.5 

Minimum  outer  radius  (cm) 

125.5 

124.5 

124.5 

Field  ARf  (G) 

-590.0 

-590 

-5900 

age  fields  have  been  ignored  and  therefore  the  ring 
current  is  not  a  relevant  parameter.  Also  at  this  high  y 
the  beam  minor  diameter  is  expected  to  be  only  a  few 
millimeters. 

In  run  267,  the  pulsed  agitator  was  turned  on  after  the 
ring’s  major  radius  became  121  cm.  Figure  3(a)  shows 
the  radial  excursion  of  a  typical  electron  that  was  located 
at  at  the  turning  on  of  the  agitator.  After  a  single 
pass  through  the  agitator  the  electron  obtains  enough  ra* 
dial  excursion  to  enter  the  extractor  and  is  extracted. 
Figure  3(b)  shows  that  the  elearon  at  the  disturbance 
obtains  a  transverse  velocity  approximately  2.8xlO~V 
Equation  (2)  predicts  a  Aur“‘2.7x  10“V  In  addition, 
the  numerical  results  show  that  the  electron  gyrates 
around  Bt  with  a  I -cm  radius,  which  is  also  the  radius 
predicted  by  Eq.  (3). 

In  run  266,  the  electron  started  at  1 10  cm  and  was 
moved  radially  by  the  mismatching  field.  The  elapsed 
time  from  the  minor  axis  to  the  agitator  is  =4.3  //sec, 
that  corresponds  to  an  average  radial  velocity  of 
=  2.2x  10*  cm/s*c.  The  amplitude  of  the  slow  mode  is 
less  than  2  mm.  Figure  4(a)  shows  the  radial  excursions 
of  a  typical  electron  in  the  r,0  plane  and  Fig.  4(b)  shows 
a  top  view  of  its  orbit.  The  electrons  reach  the  extractor 
with  a  vertical  displacement  from  the  midplane  that  is 
only  a  few  millimeters.  For  the  reason  given  in  the  pre¬ 
vious  section,  the  run  was  terminated  after  the  electron 
propagated  =  30  cm  inside  the  extractor. 

In  run  268,  y  was  increased  to  4(X)  with  a  correspond¬ 
ing  increase  in  the  value  of  magnetic  fields.  Figure  S(a) 


shows  the  radial  excursions  of  the  electron  and  Fig.  5(b) 
is  a  top  view  of  the  orbit.  The  coherence  of  the  radial 
excursions  is  remarkable.  We  have  found  that  this 
coherence  is  preserved  even  when  Eq.  (1)  is  not  satisfied 


0  0.4  os  1.2  1.6  2.0 


B/yr 


FIG.  3.  (a)  Radial  excursions  of  a  typical  electron  and  (b) 
its  corresponding  normalized  transverse  velocity  for  the  run 
267. 
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FIG.  4.  (a)  Radial  excursions  of  a  typical  electron  and  <b) 
top  view  of  its  trajectory  in  the  x,y  plane  for  run  266. 


exactly,  i.c.,  when  the  value  of  field  is  off  by  tens  of 
gauss.  Finally,  by  operation  at  /“I  or  2  instead  of  at 
/”3,  the  value  of  Bt  can  be  substantially  reduced. 

Several  successive  runs  with  y“400  electrons  that  arc 
located  at  different  positions  on  the  minor  cross  section 
of  a  3-mm-diam  ring  show  that  all  these  electrons  arrive 
at  the  entrance  of  the  extractor  with  less  than  1  nsec 
time  difference.  In  addition,  the  results  show  that  al¬ 
though  the  vertical  width  of  the  ring  has  been  slightly  re¬ 
duced,  its  corresponding  radial  width  has  been  increased 
by  approximately  a  factor  of  3.  This  spread  of  the  beam 
is  probably  associated  with  the  slightly  different  value  of 
AB:  experienced  by  the  different  electrons  as  a  result  of 
the  finite  gradient  of  the  agitator. 

In  conclusion,  we  have  developed  a  new  extraction 
scheme  that  is  practical  and  has  the  potential,  since  all 
the  electrons  of  the  ring  perform  coherent  motion,  to 
lead  to  a  very  high  extraction  efficiency. 

The  authors  are  grateful  to  Professor  D.  Kerst  and  Dr. 
J.  Golden,  Dr.  P.  Sprangle,  and  Dr.  R.  Faehl  for  many 
illuminating  discussions.  This  work  was  supported  by 
the  Space  and  Naval  Warfare  Systems  Command  and 
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Flo.  5  (a)  Radial  excursion  of  a  typical  electron  and  (b)  top 
view  of  its  trajectory  in  the  x,y  plane  for  run  268. 
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Improved  Beam  Confinement  in  the  Modified  Betatron  with  Strong  Focusing 

C.  A.  Kapetanakos,  L.  K.  Len,^*’  T.  Smith,  J.  Golden,***’  K.  Smith,*'’  S.  J.  Marsh,*'’  D.  Dialetis,**” 

J.  Mathew,  P.  Loschialpo,  and  J.  H.  Chang 
Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  D.  C.  20S75 
(Received  16  January  1990) 

The  lifetime  of  the  circulating  electron  beam  in  the  Naval  Research  Laboratory  modified  betatron  has 
been  increased  by  more  than  2  orders  of  magnitude  with  the  addition  of  strong  focusing  windings  to  the 
device.  The  injected  beam  consistently  spirals  from  the  injector  to  the  minor  axis  and  is  trapped.  The 
— O.S-kA  trapped  electron  ring  has  been  accelerated  above  10  MeV  from  the  injr^ion  energy  of  O.S  to 
0.6  MeV.  The  beam  acceleration  has  been  confirmed  not  only  by  the  x-ray  attenuation  technique  but 
also  with  the  detection  of  photoneutrons. 

PACS  numbers:  4l.80.Ec.  29.20.Fj 


Accelerators  that  combine  high-current  capability  and 
high  effective  accelerating  gradient  arc  currently  under 
development  in  several  laboratories.'*’  These  accelera¬ 
tors  have  two  common  features,  namely,  strong  focusing 
that  improves  their  current-carrying  capability  and  recir¬ 
culation  that  enhances  their  eff'ective  accelerating  gra¬ 
dient. 

Among  these  accelerators  is  the  modified  betatron. 
This  device  is  currently  under  investigation  at  the  Uni¬ 
versity  of  California,  Irvine^  and  also  at  the  Naval 
Research  Laboratory'  (NRL).  fn  its  original  form  the 
modified  betatron  consists  of  a  strong  toroidal  magnetic 
field  and  a  time-varying  betatron  field  that  is  responsible 
for  the  acceleration. 

Extensive  studies  of  beam  capture  and  confinement  in 
the  NRL  modified  betatron  led  to  the  formation  of  elec¬ 
tron  rings  with  circulating  current  as  high  as  3  kA  and 
have  furnished  valuable  information  on  the  critical  phys¬ 
ics  issues  of  the  concept.  In  addition,  these  studies  have 
shown  that  over  a  wide  range  of  parameters  the  ring  life¬ 
time  was  limited  to  a  few  microseconds  which  is  compa¬ 
rable  to  the  magnetic-field  diffusion  time  through  the 
vacuum  chamber.  Thus,  it  became  apparent  from  these 
results  that  the  magnetic-field  configuration  of  the  device 
had  to  be  modified  to  increase  the  beam  lifetime  and 
thus  to  achieve  acceleration.  A  decision  was  made  to 
proceed  rapidly  with  the  design,  fabrication,  and  installa¬ 
tion  of  strong  focusing  windings. 

Stellarator’  (four  twisted  windings  that  carry  current 
in  alternate  directions)  and  Torsatron'®  (two  twisted 
windings  that  carry  current  in  the  same  direction)  wind¬ 
ing  configurations  were  considered.  Both  configurations 
have  advantages  and  shortcomings.  The  stellarator  con¬ 
figuration  was  finally  selected  not  only  oecause  of  the 
small  net  vertical  field  and  the  lower  current  per  winding 
but  also  because  it  is  compatible  with  our  contemplated 
extraction  scheme. " ■ '  ^ 

This  paper  briefly  describes  the  NRL  modified  beta¬ 
tron  in  its  latest  form  and  summarizes  the  most  impor¬ 


tant  experimental  results  after  the  installation  of  the 
strong  focusing  windings.  Although  experiments  were 
carried  out  at  various  background  pressures,  most  of  the 
results  reported  in  this  paper  were  limited  to  pressures 
between  2xl0~*  and  8x10”®  Torr.  In  this  pressure 
range  the  electron-beam  electrical  neutralization  time  by 
direct  coliisional  ionization  of  the  background  gas  has 
been  estimated  to  be  between  180  and  45  psec.  Work  is 
in  progress  to  reduce  the  background  pressure  by  at  least 
an  order  of  magnitude  and  thus  avoid  substantial  plasma 
formation  over  the  entire  beam  lifetime.  This  is  in  con¬ 
trast  to  the  main  effort  in  the  Stellatron  at  Irvine^  that  is 
focused  on  the  formation  of  runaway  electron  beams  and 
therefore  the  pressure  is  incidentally  high. 

Description  of  the  experiment. —  (1)  External  magnet¬ 
ic  fields:  The  NRL  modified  betatron  comprises  three 
different  externa!  magnetic  fields:  the  betatron  field  that 
is  a  function  of  time  and  is  responsible  for  the  accelera¬ 
tion  of  the  electrons,  the  toroidal  magnetic  field  that 
varies  only  slightly  during  the  acceleration  of  the  elec¬ 
tron  ring,  and  the  strong  focusing  field  that  also  has  a 
very  weak  time  dependence.  The  coils  that  generate 
these  three  fields  are  shown  schematically  in  Fig.  1. 

(2)  The  betatron  field;  The  betatron'’  magnetic  field 
controls  mainly  the  major  radius  of  the  gyrating  electron 
ring  and  is  produced  by  eighteen  air-core,  circular  coils 
connected  in  series.  Their  total  inductance  is  approxi¬ 
mately  530  pH.  The  coils  arc  powered  by  an  8.64-mF 
capacitor  bank  that  can  be  charged  up  to  17  kV,  At  full 
charge,  the  bank  delivers  to  the  coils  a  peak  current  of 
about  45  kA.  The  current  flowing  through  the  coils  pro¬ 
duces  a  field  that  varies  sinusoidally  with  time  having  a 
quarter  period  of  2.5  msec  and  an  amplitude  on  the 
minor  axis  at  peak  charging  voltage  equal  to  2.1  kG. 
Immediately  after  the  peak  the  field  is  crowbarred  with  a 
4.5-msec  decay  time. 

The  flux  condition  and  field  index  are  adjusted  by  two 
sets  of  trimmer  coils  that  arc  connected  in  parallel  to  the 
main  coils.  The  current  through  the  trimmers  is  adjust- 


2374 


Work  of  the  U.  S.  Government 
Not  subject  to  U.  S.  copyright 


009 


Volume  64,  Number  20 


PHYSICAL  REVIEW  LETTERS 


14  May  1990 


Centrot 


FIG.  1.  Elevation  of  the  device  showing  the  vertical  field 
(V.F.),  toroidal  field  (T.F.),  stellarator  windings,  and  some 
structural  components. 


cd  with  series  inductors.  Typically  ~107c-l5%  of  the 
total  current  flows  through  the  trimmers. 

(3)  The  toroidal  field:  The  toroidal  magnetic  field 
controls  mainly  the  minor  cross  section  of  the  electron 
ring  and  the  growth  rate  of  several  unstable  collective 
modes.  This  field  is  generated  by  twelve  air-core,  rec¬ 
tangular  coils  that  are  connected  in  series. 

The  coils  are  made  of  aluminum  square  tubing  and 
have  a  150  cm  height  and  135  cm  width.  The  total  in¬ 
ductance  of  the  twelve  coils  is  —85  uH  and  are  powered 
by  a  24-mF  capacitor  bank  that  can  be  charged  to  a 
peak  voltage  of  10.6  kV.  At  peak  voltage,  the  bank 
delivers  to  the  coils  —180  kA.  This  current  produces  a 
field  that  varies  sinusoidally  with  time  having  a  quarter 
period  of  1.9  msec  and  an  amplitude  on  the  minor  axis 
equal  to  3.9  kG. 

(4)  The  strong  focusing  field:  The  strong  focusing 
field  improves  the  confining  properties  of  the  other  two 
fields  by  reducing  the  sensitivity  of  the  beam  centroid 
and  individual  electrons  to  energy  mismatch  and  energy 
spread. 

In  the  NRL  modified  betatron  the  strong  focusing 
field  is  generated  by  four  twisted  windings  carrying 
current  in  alternate  directions.  The  left-handed  windings 
arc  located  23.4  cm  from  the  minor  axis  and  have  a 
209.4-cm  period,  i.e.,  there  arc  three  periods  over  the  cir¬ 
cumference  of  the  torus.  They  are  supported  by  epoxy- 
reinforced  graphite  jackets  and  have  been  designed  to 
carry  up  to  25  kA  The  windings  arc  connected  in  scries 
and  the  current  temporal  profile  is  controlled  by  a  ballast 
inductor 

The  stellarator  field  index  is  defined,  in  the  linear 


approximation,  as'° 

ny  —  —  (nJ*e„)^ar(i/2n,oft«o . 

where 

—  ^odopo^ !i(2opo) , 

no"  |e  |Uo/m,  a— 2a/L,  tq  is  the  major 

radius,  AfJ  is  the  derivative  of  the  modified  Bessel  func¬ 
tion,  and  n^o  and  n«o  are  the  cyclotron  frequencies  of 
the  vertical  and  toroidal  fields  on  the  minor  axis.  Since 
in  the  NRL  modified  betatron  /«,  n>o.  and  H-o/y 
remain  approximately  constant  during  acceleration, 
scales  inversely  proportional  to  the  relativistic  factor  y. 
The  stellarator  field  index  is  related  to  the  maximum 
gradient  g  I  "(dR;/dr of  the  stcl- 
laralor  field  by  n„  “  —  g^rfJlaBsoBgQ.  For  “2.0  kG, 
5jo“25  G.  o“— 0.03/cm,  po “23.42  cm.  and  /„“I9 
k A,  rin  ==  1 4,  and  g  —  20  G/cm. 

(5)  The  vacuum  chamber  and  diode:  The  !00-cm  ma¬ 
jor  radius,  15.2-cm-insidc  minor  radius  vacuum  chamber 
has  been  constructed  using  epoxy-reinforced  carbon 
fibers.  The  desired  conductivity  is  obtained  by  embed¬ 
ding  a  phosphor  bronze  screen  inside  the  body  of  the 
graphite.  The  graphite  is  2.5  mm  thick  and  has  a  sur¬ 
face  resistivity  of  26.6  mn  on  a  square.  The  screen  has 
250x250  wires  per  inch  and  is  ';adc  of  40-pm-diam 
wire  with  an  equivalent  surface  resistivity  of  12.8  mfl  on 
a  squ&TC.  The  calculated  resistance  for  the  entire  vacu¬ 
um  chamber  is  57  mfl.  The  measured  dc  resistance  of 
the  toroidal  vacuum  vessel  is  68 ±2  mfl.  The  outside 
surface  of  the  chamber  is  covered  with  a  6.3-mm-thick. 
epoxy-reinforced  fiberglass  layer. 

This  novel  construction  technique  has  several  attrac¬ 
tive  features,  including  controllable  resistivity  and  thus 
magnetic-field  penetration  time,  high  stiffness,  and  ten¬ 
sile  strength,  high-radiation  resistance  (up  to  5(X)  Mrad) 
and  low-outgassing  rate  (  — 10~*  Torr/secem^). 

The  electrons  arc  emitted  from  a  thin  carbon  disk 
mounted  at  one  end  of  a  2.5-cm-diam  cylindrical  cathode 
stalk  clad  with  molybdenum.  The  carbon  disk  is  flush 
with  the  molybdenum  cladding  and  faces  the  circular 
opening  of  the  conical  anode,  that  is  located  8.7  cm  from 
the  minor  axis  of  the  toroidal  chamber.  To  minimize  the 
magnetic-field  component  transverse  to  the  emitting  sur¬ 
face  of  the  cathode  the  core  of  the  cathode  stalk  is  made 
out  of  high  ap  ferromagnetic  material.  The  unsaturated 
state  of  the  ferromagnetic  material  is  prolonged  by  a  thin 
copper  housing. 

Experimental  results.  —  During  the  first  microsecond 
following  injection,  the  beam  transverse  motion  is  stud¬ 
ied  by  monitoring  the  light  emitted  from  a  10-pm-thick 
polycarbonate  foil  that  is  stretched  across  the  minor 
cross  section  of  the  vacuum  chamber. '  The  foil  is  graph¬ 
ite  coated  on  the  upstream  side  to  avoid  electrostatic 
charging.  Figure  2  shows  two  open-shutter  photographs 
of  the  light  emitted  as  the  ring  passes  through  the  foil 
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FIG.  3.  Output  from  a  Rogowski  coil  that  is  located  15° 
downstream  of  the  diode.  The  plateau  (trapped  current)  corre¬ 
sponds  to  525  A. 
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FIG.  2.  Open-shutter  photographs  of  the  light  emitted  when 
the  beam  passes  a  IO->im-thick  foil. 


The  light  spots  near  the  edges  of  the  photograph  are  po¬ 
sition  fiducials  produced  by  a  ring  of  light-emitting 
diodes  located  approximately  0.7  cm  from  the  vacuum 
chamber  wall.  Reflections  from  the  wall  are  also  notice¬ 
able  near  the  upper  edge  of  the  photograph. 

For  several  combinations  of  injection  parameters  the 
beam  consistently  spirals  from  the  injection  position  to 
the  minor  axis  and  is  trapped.  The  transit  time  of  the 
beam  to  the  minor  axis  is  typically  less  than  0.5  //sec.  In 
shot  No.  5172  both  the  beam  current  and  Be  are  low  and 
thus  the  bounce  frequency  is  high.  As  the  beam  spirals 
from  the  diode  to  the  minor  axis  it  creates  distinct  light 
spots  as  it  passes  through  the  foil.  By  increasing  the 
beam  current  and  Be  the  bounce  frequency  is  reduced 
and  thus  the  light  spots  start  to  overlap  and  the  spiral  be¬ 
comes  continuous  as  in  shot  No.  5214. 

The  trapping  mechanism  is  presently  unknown.  A 
likely  candidate  is  the  wall  resistivity.  However,  the  pre¬ 
dicted  decay  rate  r~'  from  the  linear  theory'^  for  the 
parameters  of  the  experiment  is  between  10  and  20  psec, 
i.e  ,  too  long  to  explain  the  experimental  results 

The  circulating  electron  ring  current  is  monitored  with 
two  Rogowski  coils  that  are  located  inside  the  vacuum 
chamber.  A  typical  wave  form  is  shown  in  Fig.  3.  Al¬ 
though  some  losses  occur  immediately  after  injection,  the 
circulating  current  remains  practically  constant  during 
the  first  10  /rsec.  Observations  are  limited  to  10  ^sec  by 
the  signal  level  and  the  100-^/sec  passive  integrator.  At¬ 
tempts  to  use  active  integrators  have  been,  so  far,  unsuc¬ 
cessful. 

In  the  results  reported  here  the  injected  electron  beam 
is  produced  by  a  diode  with  a  l-cm-diam  carbon  cathode 


disk  that  is  matched  to  a  l-cm-diam  anode  hole  For 
such  a  beam  radius  the  maximum  trapped  current  is  be¬ 
tween  0.5  and  0.6  kA.  The  results  indicate  that  the 
trapped  current  increases  as  the  beam  diameter  in¬ 
creases.  Experiments  with  larger  radius  beams  are 
presently  in  progress. 

The  ring  lifetime  is  inferred  from  the  x  rays  produced 
when  the  beam  strikes  a  2.5-cm-widc,  1-mm-thick  lead 
limiter.  The  x  rays  are  monitored  by  a  collimated  x-ray 
detector  (scintillator-photomultiplier  tube)  that  is 
housed  inside  a  lead  box.  In  the  results  reported  in  this 
paper,  the  x  rays  enter  the  scintillaior  through  a  1.94- 
cm-diam  hole  and  the  detector  is  located  10.8  m  from 
the  target.  A  typical  wave  form  of  the  x-ray  monitor  is 
shown  in  Fig.  4.  From  the  value  of  the  B-  field  it  can  be 
easily  computed  that  the  main  peak  of  the  x-ray  signal 
corresponds  to  a  particle  energy  near  1 1.0  MeV.  The  en¬ 
ergy  of  the  electrons  that  are  lost  at  a  later  time  is  obvi¬ 
ously  higher. 

The  x-ray  signal  is  very  reproducible  in  both  ampli¬ 
tude  and  shape.  For  the  first  200-300  /rsec  wc  do  not 
observe  any  x  rays.  The  time  at  which  x  rays  arc  initial¬ 
ly  observed  and  the  times  the  various  signal  peaks  occur 
are  directly  proportional  to  Be  and  inversely  proportional 
to  dB-Jdt.  In  addition,  the  temporal  occurrence  of  the 
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x-ray  signal  appears  to  be  independent  of  the  trapped 
current  for  at  least  up  to  0.5  kA  and  also  of  the  back¬ 
ground  pressure  up  to  6 x  10 Torr.  However,  the  am¬ 
plitude  of  the  x-ray  signal  decreases  rapidly  with  increas¬ 
ing  pressure  above  8x10”*  Torr.  For  the  results  of  Fig. 
4,  the  electrons  gain  energy  at  the  rate  of  0.39  kcV/ 
revolution.  Those  electrons  lost  at  the  main  x-ray  peak 
have  performed  in  excess  of  28000  revolutions  around 
the  major  axis.  In  our  best  results  the  corresponding 
number  is  ~ 35000  revolutions  and  the  electron  energy 
about  1 2  MeV. 

The  spiky  shape  of  the  x-ray  signal  and  therefore  the 
beam  loss  is  consistent  with  the  excitation  of  the  cyclo¬ 
tron  resonance.'*  '*  Figure  4  also  shows  the  ratio  BJB, 
versus  time.  The  solid  circles  are  from  the  experiment 
and  the  crosses  are  the  prediction  of  the  theory.  The 
poor  argument  between  theory  and  experiment  for  low-/ 
values  is  probably  associated  with  the  omission  of 
toroidal  effects  in  the  theoretical  model.  The  results  are 
also  consistent  with  the  excitation  of  the  electron-cyclo¬ 
tron  instability.'^  This  instability  is  caused  by  the  cou¬ 
pling  of  the  electron-cyclotron  mode  to  the  TEu 
waveguide  mode  of  the  torus.  Work  is  in  progress  to  pin¬ 
point  the  exact  cause  of  the  beam  loss. 

The  x-ray  signal  has  been  studied  with  the  usual  at¬ 
tenuation  technique.  A  layer  of  lead  more  than  5  cm 
thick  is  needed  to  completely  eliminate  the  signal.  A 
1.3-cm-thick  layer  does  not  have  any  effect  on  the  signal 
except  on  the  initial  spike  that  is  due  to  the  injected 
beam.  In  addition  to  the  x-ray  attenuation  technique, 
the  beam  acceleration  has  been  confirmed  with  the  pho¬ 
toneutrons  produced  from  the  D(y,n)H  reaction.  A 
plastic  tube  in  the  form  of  a  ring  containing  heavy  water 
was  inserted  behind  the  limiter.  The  photons  produced 
on  the  target  photodisintegrate  the  deuterium  and  pro¬ 
duce  neutrons  that  are  monitored  with  a  rhodium  activa¬ 
tion  detector.'*  The  number  of  counts  measured  during 
the  first  minute  exceeds  7  standard  deviations. 

The  intensity  and  polarization  of  the  synchrotron  radi¬ 
ation  has  been  computed  numerically.  At  low  energy 
(<  1  MeV)  the  radiation  spectrum  is  dominated  by  a 
peak  at  the  Doppler-shifted  B*  cyclotron  frequency.  As 
the  electron  energy  increases  the  effect  of  the  cyclo¬ 
tron  motion  is  reduced  and  the  spectrum  approaches  that 
of  a  purely  circular  orbit. 

Two  high-gain  (50  db)  heterodyne  receivers  are  used 
to  monitor  the  emitted  radiation  in  the  bands  26.5-28.5 
and  36.5-38.5  GHz  during  acceleration.  The  measure¬ 
ments  show  that  the  amplitude  of  the  radiation  scales 
linearly  with  the  trapped  ring  current.  During  the  first 
200  nsec,  both  the  amplitude  and  polarization  are  con¬ 
sistent  with  the  predicted  single-particle  emission.  In  ad¬ 
dition,  the  temporal  shift  of  the  first  peak  is  in  good 
agreement  with  the  predicted  toroidal  field  and  energy 


dependence  of  the  Doppler-shifted  cyclotron  frequency. 
However,  the  late-time  behavior  of  the  radiation  signals 
is  not  presently  well  understood.  The  radiation  signal 
starts  to  decay  after  ~  200  pscc,  contrary  to  the  predic¬ 
tion  of  the  theory.  It  is  likely  that  electron  losses,  plasma 
formation,  or  beam  displacement  out  of  the  field  of  view 
of  the  detectors  are  responsible  for  the  observed  dis¬ 
crepancy. 

Tbe  authors  arc  grateful  to  Professor  D.  Kerst,  Dr.  P. 
Sprangle,  and  Dr.  Y.  H.  Seo  for  many  illuminating  dis¬ 
cussions.  This  work  was  supported  by  the  ONR  and  the 
Space  and  Naval  Warfare  Systems  Command. 

Note  added. — Since  the  submission  of  the  paper,  the 
trapped  current  increased  to  0.9- 1.0  kA  and  the  beam 
energy  to  15-16  MeV. 
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Compact,  high-current  accelerators  and  their  prospective  applications* 

C.  A.  Kapetanakos/  L.  K.  Len,*’  T.  Smith.  D.  Dialetis.*”  S,  J.  Marsh, P.  Loschiaipo, 
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This  paper  briefly  surveys  the  three  compact,  high-current  accelerators  that  are  presently 
under  development  in  the  United  States  in  support  of  a  national  program.  In  addition,  it 
reports  recent  experimental  results  from  the  Naval  Research  Laboratory  (NRL)  modified 
betatron  (Phys.  Rev.  Lett.  64,  2374  ( 1990)]  with  emphasis  on  the  electron-cyclotron 
resonance  that  presently  limits  the  energy  of  the  beam  to  approximately  18  MeV.  Finally,  it 
briefly  addresses  selective  existing  and  prospective  applications  of  accelerators. 


I.  INTRODUCTION 

Over  the  last  60  years  the  technology  of  high-energy 
accelerators  has  advanced  to  a  remarkable  level  of  sophisti¬ 
cation.'  These  conventional  accelerators  are  designed  to 
operate  reliably  at  high  energies  but  at  relatively  low  current, 
primarily  to  avoid  complications  related  to  beam  self-field 
effects.  Thus  the  beam  dynamics  in  conventional  accelera¬ 
tors  is  determined  mainly  by  the  externally  applied  fields. 

Over  the  last  several  years,  it  has  become  apparent  that 
accelerators  with  both  high  current  and  high-average  power 
could  have  interesting  applications  not  only  in  national  de¬ 
fense^  but  also  in  several  areas  of  the  civilian  economy.' 

In  contrast  with  the  conventional  accelerators  that 
operate  in  the  single-particle  regime,  high-current  accelera¬ 
tors  operate  in  an  uncharted  territory,  in  which  space  charge 
effects  from  the  self  and  image  fields  are  as  important  as 
externally  applied  fields.  The  free  energy  stored  in  the  in¬ 
tense  self  fields  of  these  beams  can  drive  collective  instabili¬ 
ties  and  strong  magnetic  fields  are  needed  to  reduce  the  rapid 
growth  rate  of  these  unstable  modes.  Several  of  the  key  phys¬ 
ics  issues  associated  with  the  high-current  accelerators  are 
similar  to  those  in  plasma  physics  Therefore  it  is  not  acci¬ 
dental  that  the  majority  of  researchers  that  are  presently 
developing  this  technology  have  been  trained  in  plasma 
physics  and  not  in  conventional  accelerator  technology. 

Compact,  high-current  accelerators  are  currently  under 
development  in  several  laboratories.  '  ^  As  a  result  of  their 
compactness  and  light  weight,  these  devices  are  expected  to 
be  substantially  less  costly  than  conventional  accelerators  of 
comparable  power.  In  addition,  because  of  their  high-cur¬ 
rent  carrying  capability,  these  accelerators  have  the  poten¬ 
tial  to  generate  very  powerful  electron  beams  that  can  pro¬ 
vide  high  dose  rates  at  reduced  unit  irradiation  cost. 

During  the  last  several  years  several  applications  have 
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been  suggested,'  *  such  as  water  purification,  sludge  disin¬ 
fection.  cracking  of  crude  oil,  and  excavation  of  tunnels,  all 
requiring  high-average  power  and  low  unit  irradiation  cost 
These  applications  are  presently  on  hold  waiting  for  the  de¬ 
velopment  of  inexpensive  beams. 

In  this  paper,  we  briefly  survey  the  three  compact  accel¬ 
erators  that  are  currently  under  development  in  support  of 
the  DARPA/services  program  and  report  recent  experi¬ 
mental  results  from  the  NRL  modified  betatron  The  last 
section  briefly  addresses  existing  and  prospective  applica¬ 
tions  of  accelerators. 

II.  REVIEW  OF  COMPACT,  HIGH-CURRENT 
ACCELERATORS 

Compact,  high-current  accelerators  are  currenMy  under 
development  at  the  University  of  California,  Irvine,'*  Uni¬ 
versity  of  New  Mexico.'  Pulse  Sciences  Inc.*  (PSD.Sandia 
National  Laboratory’  (SNL),  and  the  Naval  Research  Lab¬ 
oratory'  (NRL).  All  these  concepts  are  induction  accelera¬ 
tors.  The  reason  is  that  induction  accelerators  are  inherently 
low-impedance  devices  and  thus  are  ideally  suited  to  drive 
high-current  beams.  The  acceleration  process  is  based  on  ihe 
electric  field  produced  by  a  lime-varying  magnetic  field.  The 
electric  field  can  be  either  continuous  or  localized  along  the 
accelerating  path. 

This  section  briefly  reviews  the  three  electron  accelera¬ 
tion  concepts  that  are  part  of  the  national  program:  The  spi¬ 
ral  line  induction  accelerator  (SLlA)  developed  at  PSI.  the 
recirculating  linear  accelerator  (RLA)  developed  at  SNL. 
and  the  modified  betatron  accelerator  (MBA)  developed  at 
NRL.  In  addition  to  being  inductive,  these  three  accelerators 
utilize  strong  focusing  m  order  to  improve  their  current  car¬ 
rying  capability  and  recirculation  to  enhance  their  accelerat¬ 
ing  gradient. 

A.  Spiral  line  induction  accelerator* 

The  SLI A  concept  is  shown  schematically  m  Fig  !  ( a ) 
The  beam  is  injected  at  one  end  of  the  spiral  hue  and  is  accel¬ 
erated  at  the  advanced  test  accelerator  ( .ATA  )  type  acceler¬ 
ating  cells.  The  minor  radius  of  the  beam  is  controlled  by  an 
axial  magnetic  field  and  the  beam  is  guided  at  the  bends  by  a 
static  vertical  and  strong  focusing  field  At  the  exit,  the  beam 
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opment  program 


energy  is  equal  to  the  injection  energy  plus  the  energy  gain 
per  pass  times  the  number  of  passes. 

Since  the  SLIA  is  open-ended,  difficulties  associated 
with  injection  and  extraction  are  avoided.  Funhermore,  this 
configuration  provides  flexibility  in  the  selection  of  shape  of 
the  electron  pulses.  Important  outstanding  technical  issues 
associated  with  the  SLIA  concept  are  the  excitation  of  the 
beam  breakup  instability  (BBU),  the  matching  of  fields  be¬ 
tween  the  straight  sections  and  the  bends,  and  the  develop¬ 
ment  of  high-frequency  power  supply  for  driving  and  reset¬ 
ting  the  ferromagnetic  cores. 

A  proof-of-concept  (POC)  experiment  is  presently  un¬ 
der  design  and  construction.  Its  most  important  parameters 
are  listed  in  Table  I.  The  3.5  MeV,  10  kA,  35  nsec  duration 
injected  beam  pulse  will  be  accelerated  to  9.5  MeV  by  pass¬ 
ing  twice  through  the  two,  1.5  MeV  accelerating  units 

B.  Recirculating  linear  accelerator^ 

The  RLA  is  shown  schematically  in  Fig.  1  (b).  The  in¬ 
jected  beam  is  accelerated  by  dielectric  cavities  that  are 
placed  in  tandem  along  the  straight  sections  of  the  race 
track.  Beam  focusing  is  provided  by  an  ion  channel  formed 
by  a  low-energy,  low-current  electron  beam.  A  time-rising 
vertical  magnetic  field  guides  the  beam  at  the  bends. 

Important  advantages  of  RLA  concept  are  the  suppres¬ 
sion  of  BBU  and  the  wide  energy  bandwidth  provided  by  the 
ion  channel.  In  addition  to  injection  and  extraction,  main¬ 
taining  the  integrity  and  centering  of  the  ion  channel  during 
acceleration  and  avoiding  beam  erosion  are  key  technical 
issues. 

A  POC  experiment  is  currently  under  construction.  Its 
most  important  parameters  are  listed  in  Table  I.  The  3.5 
MeV,  10  kA,  40  nsec  duration  electron  beam  pulse  will  be 
accelerated  to  8-10  MeV  by  passing  three  times  through  the 
two,  1  MeV  dielectric  cavities. 


C.  The  modified  betatron  accelerator*-” 

In  both  the  SLIA  and  the  RLA,  the  accelerating  electric 
field  is  localized  at  the  gaps.  In  contrast,  the  electnc  field  in 
the  MBA  is  continuous  along  the  entire  beam  path  and  thus 


TABLE  I  Parameter?,  of  proof-of-concept  (  POC)  experiments 


SLIA 

RLA 

MB.A 

Injected  beam 

Injected  beam 

Injected  beam 

Energy  -  3  5  MeV 

Energy  3  5  MeV 

Energy  0  b-0  7  MeV 

Current  ~  !0  kA 

Current  10  kA 

Current  ^  2  kA 

Duration  -  35  nsec 

Duration  4()  nsec 

Duration  —40  nsec 

Axia!  mag  field  5  5  kG 

Torus  major  radius  100  cm 

Vert  mag  field  ,.0  32  kG 

Channel  radius  I  5  cm 

Torus  mmor  radius  1 5  cm 

S  F  gradient  sO.25  kG/cm 

Channel  dcnsiix.-ljeam  density  0  4 

Toroidal  mag  field  2-5  kG 

No  of  passes  2 

No  of  passes  3 

Vertical  mag  field  ff.  0-2  kG 

No  of  acc,  units  2 

No  of  casities  2 

Current  in  S  F  windings  20-  25  k  A 

Energy  gain/acc  unit  1  5  MeV 

Energy  gain/caviiy  1  MeV 

No  of  S  E'  field  periods  b 

1  5  G//isec 
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the  accelerator  is  constrained  to  circular  shape.  An  artist’s 
conception  of  the  MBA  is  shown  in  Fig.  1(c).  The  injected 
beam  is  accelerated  by  the  time- varying  vertical  (betatron) 
magnetic  field  S,.  The  toroidal  magnetic  field  controls 
mainly  the  minor  radius  of  the  beam  and  the  growth  rate  of 
the  various  unstable  modes  while  the  vertical  and  strong  fo¬ 
cusing  fields  control  mainly  the  major  radius  of  the  electron 
ring.  The  stability  properties  of  the  stellarator  windings  for 
high-current  beams  were  studied  initially  by  Gluckstem”  in 
linear  geometry  and  by  Roberson'^  et  at.  in  toroidal  geome- 
try. 

Important  advantages  of  the  MBA  are  the  natural  syn¬ 
chronization  of  B.  and  particle  energy,  the  absence  of  the 
BBU  mode  and  its  compact  size  and  light  weight.  Outstand¬ 
ing  technical  issues  are  the  suppression  of  the  cyclotron  re¬ 
sonances,  collective  instabilities,  and  the  extraction  of  the 
beam. 

A  POC  experiment  is  currently  in  operation  at  NRL. 
The  various  parameters  of  the  device  are  listed  in  Table  I.  Its 
objective  is  to  study  the  critical  physics  issues  associated 
with  the  concept  and  to  accelerate  at  1  kA,  20  nsec  electron 
beam  pulse  to  20  MeV  with  subsequent  extraction  of  the 
beam.'*  Presently,  the  trapped  electron  beam  has  been  ac¬ 
celerated  to  between  17-18  MeV  from  the  injected  energy  of 
0. 6-0.7  MeV.  The  maximum  trapped  current  is  up  to  1.5  kA 
and  the  trapping  efficiency  is  as  high  as  75%. 

III.  RECENT  EXPERIMENTAL  RESULTS  FROM  THE  MBA 

The  NRL  modified  betatron  accelerator  has  been  de¬ 
scribed  previously'  and  a  summary  of  the  initial  experi¬ 
mental  results  has  been  published.''*’  In  this  paper,  we  sum¬ 
marize  the  highlights  of  the  most  recent  results. 

During  the  first  microsecond  following  injection,  the 
beam  motion  is  inferred  from  the  light  emitted  when  the 
electrons  pass  through  a  thin  (3  yum)  plastic  foil  stretched 
across  the  minor  cross  section  of  the  vacuum  chamber.  The 
open  shutter  photographs  show  that  for  a  wide  range  of  ex¬ 
perimental  parameters  and  after  fine  tuning  the  fields,  the 
injected  beam  consistently  spirals  from  the  injector  to  the 
vicinity  of  the  minor  axis  and  is  trapped.  The  beam  travel 
time  from  the  injector  to  the  minor  axis  depends  on  the  val¬ 
ues  of  the  various  fields  and  is  typically  between  300-1000 
nsec. 

The  trapping  of  the  beam  is  most  likely  due  to  the  wall 
resistivity  of  the  vacuum  chamber.  The  physical  mechanism 
of  the  resistive  trapping  can  be  understood  as  follows.  Con¬ 
sider  a  pencil-like  electron  beam  inside  a  cylindrical  pipe 
with  finite  resistivity  wall.  Becau.se  of  the  resistivity  of  the 
wall,  the  electric  and  magnetic  images  do  not  coincide  and 
thus  the  centroid  of  an  off-axis  beam  experiences  a  negative 
radial  magnetic  field  component.  This  field  component 
crossed  with  the  axial  velocity  of  the  beam  produces  a  poloi- 
dal  force,  which  is  in  the  opposite  direction  to  the  poloidal 
motion  of  the  beam  In  the  absence  of  the  strong  focusing 
and  when  the  self-fields  dominate  the  external  fields  (high- 
current  regime),  the  poloidal  force  in  cvinjunction  with  the 
axial  (toroidal)  magnetic  field  drives  the  beam  to  the  wall 
(drag  instability ) .  How  ever,  in  the  presence  of  strong  focus- 
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tng,  the  direction  of  the  poloidal  motion  can  be  reversed  and 
the  beam  spirals  to  the  minor  axis. 

The  resistive  wall  trapping  mechanism  has  been  ana¬ 
lyzed  previously."*  The  linear  theory  was  done  for  a  beam 
near  the  minor  axis.  In  this  case,  the  relevant  diffusion  time 
is  T,  ^Hoo/lp,  where  a  is  the  minor  radius  of  the  torus  and  p 
its  surface  resistivity.  For  the  parameters  of  the  NRL-MBA, 
r,  ==  10  //sec.  The  linear  decay  rate  T  “  '  from  the  theory  is 
between  10-15  //sec,  as  expected.  This  time  is  at  least  a  factor 
of  30  too  long  to  explain  the  results. 

In  contrast  to  the  analysis  that  assumes  the  beam  is  near 
the  minor  axis,  the  beam  in  the  experiment  during  injection 
is  near  the  wall  and  therefore  the  relevant  diffusion  time  is 
Tj  =  po^/v^p,  where  A  is  the  thickness  of  the  wall.  The 
ratio  r,  /r^  =  rra/^  =  300.  The  fact  that  r,  is  substantially 
shorter  than  the  observed  trapping  time  Is  not  surprising 
since  the  beam  remains  near  the  wall  only  for  a  shon  period 
of  time.  Expieriments  are  presently  under  way  to  test  the 
hypothesis  of  the  resistive  trapping. 

The  X  rays  are  monitored  by  three  collimated  x-ray  de¬ 
tectors  (scintillator-photomultiplier  tube)  that  are  housed 
inside  lead  boxes.  In  the  results  shown  in  Fig.  2,  the  x  rays 
enter  the  scintillator  through  a  1.94  cm  diam  tube  and  the 
detector  is  located  10.8  m  from  the  vacuum  chamber.  As  a 
rule,  the  shape  of  the  x-ray  signal  recorded  by  all  three  detec¬ 
tors  is  spiky  and  the  peaks  always  occur  at  the  same  value  of 
Bqq/B^  {BgQ  is  the  toroidal  and  B^  the  vertical  magnetic 
fields  on  the  minor  axis)  independently  of  the  current  flow¬ 
ing  in  the  stellarator  windings. 

In  addition  to  the  x-ray  pulse.  Fig.  2  shows  the  ratio  of 
B„„/B^,  at  the  peaks  of  the  signal  The  solid  circles  are  from 
the  experiment  and  the  cros.ses  are  from  the  resonance  condi¬ 
tion  B„„/B^  =  (2/’  -  l)/2/,  /  =  1,2,...,  of  the  cyclotron 
resonance."* The  solid  line  shows  the  trend  of  the  experi¬ 
mental  results  with  time.  For  /  =  8,  9,  and  10,  experiment 
and  theory  are  in  goc>d  agreement.  However,  for  the  remain¬ 
ing  /  values  there  is  noticeable  divergence  between  the  theo¬ 
retical  predictions  and  the  experiment. 

The  theory  of  the  c>clotron  resonance  has  been  desel- 
oped  under  the  assumption  that  the  beam  is  located  on  the 
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minor  Mis.  However,  the  experimental  observations  indi¬ 
cate  that  the  electron  ring  starts  to  move  off  axis  after  200- 
300  //sec.  Recently,  the  theory  of  the  cyclotron  resonance 
has  been  extended  to  a  beam  that  is  located  on  the  midplane 
but  off  the  minor  axis.  For  such  a  beam  the  resonance  condi¬ 
tions  becomes 

=  (2/'- »)/2/.  (I) 

where  is  the  major  radius,  the  cyclotron  frequency 
of  the  toroidal  held,  y  Is  the  relativistic  factor,  and  0^  is  the 
normalized  toroidal  velocity.  Perfect  match  between  the  ex¬ 
perimental  results  and  the  revised  theory  is  obtained  if  it  is 
assumed  that  the  equilibrium  position  of  the  beam  moves 
initially  toward  the  major  axis  and  then  away  from  it  with  a 
constant  outward  radial  speed  of  0.23  mm///sec,  starting  at 
about  460  //sec. 

In  a  modified  betatron  with  strong  focusing  windings, 
there  are  four  characteristic  transverse  modes'^  6/  ^  .  In¬ 

teger  resonances  occur  when  the  frequency  of  these  modes, 
in  the  laboratory  frame,  over  the  relativistic  cyclotron  fre¬ 
quency  of  the  vertical  field  Cljo/y  on  the  minor  axis  is  an 
integer,  i.e.,  when 

where  AT  =  ±  1,  ±  2,... . 

The  cyclotron  resonance  occurs  when /a  +  +  =^. 

and,  for  the  parameters  of  the  NRL  modified  betatron  listed 
in  Table  II,  the  cyclotron  resonances  are  shown  in  Fig.  3.  It  is 
apparent  from  this  figure  that  as  the  beam  energy  increases, 
the  beam  centroid  crosses  progressively  lower-integer  num¬ 
ber  resonances. 

In  the  general  case,  the  expressions  for  the  four  charac¬ 
teristic  modes  ^  are  very  complicated.  However,  for 
modest  beam  and  strong  focusing  winding  currents,  as  those 
found  in  the  NRL  device,  the  modes  are  considerably  simpli¬ 
fied.  Under  these  conditions,  where  m  is 

the  number  of  field  periods.  In  addition,  when  >  ^>o>  Ihe 
resonance  condition  takes  the  very  simple  form  ~l, 

where  /  =  1,2,...  .  Therefore  the  cyclotron  resonance  is  due 
to  the  coupling,  caused  by  a  field  error(s),  of  the  cyclotron 
motions  associated  with  the  toroidal  and  vertical  fields. 

The  physical  mechanism  of  the  cyclotron  resonance  can 
be  easily  understood  with  the  help  of  Fig.  4.  Consider  an 
electron  rotating  in  a  uniform  vertical  field  B,  in  the  pres¬ 
ence  of  a  toroidal  field  Bg.  When  this  electron  is  at  its  equilib¬ 
rium  position,  its  orbit  is  a  circle,  centered  around  the  major 
axis,  and  the  velocity  vector  is  tangential  to  the  circular  or- 


TABLE  II  Parameters  used  in  the  computation  ofcycloiron  modes  shown 
in  Fig,  3. 


Toroidal  magnetic  field  B„ 

5  kG 

Beam  current 

I  kA 

Torus  major  radius 

100  cm 

Tonis  minor  radius 

15  cm 

Beam  minor  radius 

1  cm 

S.F.  field  periods 

6 

FIG.  3.  Beam  centroid  integer  cyclotron  resonances  for  the  parameters  list¬ 
ed  in  Table  II. 


bit.  Now,  let  us  assume  that  a  negative  vertical  field  distur¬ 
bance  AR,  is  introduced  at  a  short  segment  of  the  orbit  with 
toroidal  half-width  A(9.  As  the  electron  crosses  the  distur¬ 
bance,  it  sees  a  reduced  B,  and  tries  to  acquire  a  larger  radi¬ 
us.  Thus,  its  velocity  vector  rotates  clockwise  and  the  elec¬ 
tron  develops  a  radial  velocity  component.  As  a  result,  the 
electron  starts  to  rotate  around  Bg  and  its  orbit  is  trans¬ 
formed  from  a  circle  into  a  helix.  When  the  resonance  condi¬ 
tion  is  satisfied,  the  electron  will  return  to  the  disturbance  on 
phase  and  it  will  acquire  additional  radial  velocity  and  the 
radius  of  the  helix  will  further  increase. 

The  combined  motion  of  the  electron  around  B,  and  Bg 
is  very  similar  to  the  motion  of  the  Earth  around  the  Sun  and 
its  axis.  Fortunately,  the  time  it  takes  the  Earth  to  orbit  the 
Sun  ( 1  yr)  over  the  time  it  takes  the  Earth  to  spin  around  its 
axis  ( 1  day)  is  365.242  199,  i.e.,  is  not  an  integer.  By  select- 
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FIG  4  Physical  mechanism  of  ihecycloiron  resonance  The  disiurbance  is 
located  at  approximately  5  o'clock  at  radius  r, 
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ing  the  ratio  year/day  #  integer,  nature  has  avoided  the 
exciution  of  a  resonance  from  graviutional  disturbances. 
As  a  consequence  of  this  selection,  the  calendar  had  to  be 
reformed  three  times  (Caesar,  Augustus,  Pope  Gregory 
XIII)  over  the  last  2(X)0  years.  However,  this  is  a  rather 
small  price  to  pay. 

Since  the  initial  observation  of  the  cyclotron  resonance 
in  the  NRL  device,  we  have  uncovered  several  field  distur¬ 
bances.  By  correcting  these  field  deficiencies,  a  substantial 
increase  in  the  beam  lifetime  has  been  observed  with  subse¬ 
quent  increase  in  the  beam  energy.  One  of  the  most  impor¬ 
tant  field  errors  was  caused  by  the  ports  of  the  vacuum 
chamber.  The  ports  disturb  the  flow  pattern  of  the  current 
induced  by  the  rising  B,  and  create  magnetic  dipole-type 
field  errors.  The  measured  field  errors  from  the  ports  and 
their  Fourier  modes  are  shown  in  Figs.  5(a)  and  5(b),  re¬ 
spectively.  Figures  5(c)  and  5(e)  show  the  normalized 
transverse  velocity  components  and  5(d)  shows  the  normal¬ 
ized  toroidal  velocity  as  a  function  of  time  from  the  numeri¬ 
cal  solution  of  orbit  equations.  It  is  apparent  that  both  trans¬ 
verse  velocity  components  increase  at  the  expense  of  the 
toroidal  velocity  as  the  beam  crosses  the  /  =  10  and  9  reson¬ 
ances. 

It  has  been  suggested  that  the  beam  loss  in  the  modified 
betatron  may  be  caused  by  the  excitation  of  the  cyclotron 
instability^®  instead  of  the  cyclotron  resonance.  This  insta¬ 
bility  results  from  the  coupling  of  the  cyclotron  mode  to  the 
TE,.,  modes  of  the  toroidal  chamber  and  is  characterized  by 
the  resonance  condition  =  (x[^rl/a^  +  p')''' 

—  where  x,'  =  1.841  and  p  is  an  integer.  Figure  6 
shows  the  cavity  modes  and  the  cyclotron  mode  for 
several  values  of  The  parameter 
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FIG  6  TE,,,,  modes  of  ihe  toroidal  chamber  and  cyclotron  mode  The 
cyclotron  mode  moves  from  left  to  right  as  the  vertical  magnetic  field  B 
rises  while  the  toroidal  magnetic  field  remains  practically  constant 


R—  —  (ro/cilgf,)(da>/di)/(dilyo/(it)  is  the  normalized 
time  rate  of  the  cyclotron  mode  caused  by  the  rising  B,  field 
The  parameter  nl'"'"  is  the  exact  stellarator  field  index.'" 
Strong  interaction  is  expected  when  the  cyclotron  mode 
crosses  the  cavity  modes  near  cutoff  because  both  R  and 
n”*'"  are  small. 

To  determine  the  importance  of  the  cyclotron  instabil¬ 
ity,  an  attempt  was  made  to  measure  both  the  poloidal  B'J, 
and  toroidal  magnetic  field  components  of  the  electro¬ 
magnetic  modes  inside  the  toroidal  chamber  using  a  fast, 
broadband  B  probe.  The  results  of  these  measurements  indi¬ 
cate  that  the  peak  of  Bg  is  less  than  3x  10  '  ^  G  near  cutoff 


Field  Disturbonces  from  the  Ports 
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FIG  5  Crossing  of/ =  lOand  cyclo¬ 
tron  resonance  dunng  acceleration 
(a)  Measured  field  disturbances,  from 

the  pons  of  the  vacuum  chamber,  (hi 
Fourier  modes  of  the  disturbances.  lci 
and  (e)  normalized  transverse  velix- 
ity.  and  (d)  normalized  toroidal  velix- 
ity  versus  time 
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FIG  7.  Two  x-ray  signals  taken  under  identical  conditions  except  for  the  time  rate  of  the  vertical  field  B.  The  electron  ring  lifetime  varies  approximately 
inverse  proportionally  to  B.. 


(580  MHz),  i.e.,  the  amplitude  of  the  electromagnetic 
modes  inside  the  torus  is  negligibly  small.  Since  the  re¬ 
quired  to  excite  the  cyclotron  instability  is  between  3-10  G, 
it  is  very  unli^'-'y  thr,*  the  beam  !oes  is  due  to  the  cyclotron 
instability.  In  addition,  the  absence  of  electromagnetic 
modes  with  appreciable  amplitude  inside  the  toroidal  cavity 
precludes  the  parametric  instability*'  as  a  serious  beam  loss 
mechanism. 

The  hypothesis  that  the  beam  loss  is  due  to  the  cyclotron 
resonance  is  also  supported  by  the  fact  that  for  a  fixed  toroi¬ 
dal  magnetic  field,  the  ring  lifetime  varies  approximately  in¬ 
verse  proportionally  to  the  time  rate  of  the  vertical  magnetic 
field  dB./dt.  Figure  7  shows  the  x-ray  pulses  for  two  shots 
with  identical  parameters,  including  the  beam-trapped  cur¬ 
rent,  except  for  dB./dt.  By  reducing  dB  /dt  by  a  factor  of 
1.82,  the  beam  lifetime  increased  by  approximately  the  same 
amount. 

IV.  EXISTING  AND  POTENTIAL  APPLICATIONS  OF 
ACCELERATORS 

Table  III  lists  selective  existing  and  potential  applica¬ 
tions  of  accelerators  The  first  column  shows  the  application, 
the  second  the  dose  or  dose  rate,  and  the  third  the  desired 
average  power.  The  required  dose  varies  from  a  fraction  of  a 
Gy  for  radiation  therapy  to  -  1  MGy  for  cracking  of  crude 
oil.  i.e..  it  spans  approximateh  seven  orders  of  magnitude. 
With  the  exception  of  the  last  application  listed  in  Table  III, 
the  rest  of  them  have  already  been  discussed  in  the  scientific 
literature.' 

The  radiation-induced  cross-linkage  of  polymers  is 
probably  the  most  important  industrial  application  of  accel- 
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erators  today  and  is  used  extensively  by  several  companies 
both  in  the  U.S.  and  abroad.  For  example,  theCryovac  Divi¬ 
sion  of  W.  R.  Grace  utilizes  in  excess  of  40  electron  accelera¬ 
tors  with  an  aggregate  beam  power  in  excess  of  1  MW.  The 
annual  value  of  irradiated  products  is  —$800  M.  Raychem 


TABLE  III  Applicaiior  of  acccleraiors. 


Application 

Dose  or 
dose  rale 

Desired 

average 

power" 

( 1 )  Radiation  therapy 

0.25-1.4  Gy/min 

VL 

(2)  Sterilization 

Pbarmaceulicals 

20-30  kGy 

H 

Human  transplants 

10-20  kGy 

L 

Musical  instruments 

10  kGy 

L 

Potable  water 

1  kGy 

H 

(3)  Radiative  processing 

Sludge  disinfection 

4-10  kGy 

H 

Food  preservation 

1  0-25  kGy 

H 

Cellulose  depolymerization 

5-10  kGy 

H 

Graft  copolymerization 

10-20  kGy 

H 

Curing  of  coalings 

20-50  kGy 

M 

Vulcanization  of  silicones 

50- 150  kGy 

H 

Cross-linkage  of  jiolymers 

100-300  kGy 

H 

Vulcanization  of  rubber 

100-300  kGv 

H 

Flue-gas  clean  up 

20-30  kGv 

H 

Cracking  of  crude  oil 

1000  kGy 

H 

Radiography 

10-100  Gy/mm 

M 

Excavation  of  tunnels 

100-200  kGy 

VH 

Burning  of  nuclear  wastes 

VH 

*  VH  >  lens  of  MW,  H  =  scvtral  hundred  kW,  M  -  10-100  kW,  L  ~  kW, 
VL<kW 
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Corporation  utilizes  more  than  20  accelerators  in  the  energy 
range  0.3-10  MeV.  The  annual  value  of  the  irradiated  prod¬ 
ucts  is  ~  $1000  M.  Accelerators  are  also  used  extensively  by 
Western  Electric  to  irradiate  telephone  wires.  However,  the 
company  is  unwilling  to  release  any  information  related  to 
this  activity. 

Recently,  accelerators  have  been  utilized  to  induce  cur¬ 
ing  of  composites.^’  These  fiber-reinforced  polymers**  are 
presently  used  in  several  p-’nlications,  including  aircraft 
frames  and  wings,  automobile  oodies  and  doors,  rocket  mo¬ 
tor  cases,  oxygen  tanks,  printed  circuits,  and  tennis  rackets. 
Curing  the  composites  with  electron  beams  has  two  impor¬ 
tant  advantages;  ( 1 }  Curing  time  is  reduced  by  almost  two 
orders  of  magnitude  (from  4  days  to  2  h),  and  (2)  thermal 
stresses  are  reduced  and  thus  the  quality  of  the  product  is 
substantially  improved. 

Aerospacial  has  recently  purchased  two,  10  MeV,  20 
kW  accelerators  from  CGR  to  irradiate  rocket  motor  cases 
at  S  Mrad. 

The  last  application  listed  in  Table  III  is  presently  under 
development  at  Los  Alamos  National  Laboratory.*’  Briefly, 
it  requires  a  1 .6  GeV,  25  m A  proton  beam  that  strikes  a  Pb- 
Bi  target  producing  55  neutrons  per  proton.  The  long  half- 
life  nuclear  wastes  are  separated  from  those  with  short  half- 
life  and  are  irradiated  by  the  neutron  flux.  Transmutation 
reduces  the  half-life  of  wastes.  The  higher  atomic  weight 
actinides  form  fissionable  products  that  produce  power  in 
excess  of  that  required  to  drive  the  accelerator.  The  contem¬ 
plated  radio-frequency  quadropole  (RFQ)  linear  accelera¬ 
tor  is  1/4  mile  long,  requires  a  55  ft  wide  tunnel,  and  its  cost 
has  been  estimated  at  S500  M. 

Substantial  savings  in  the  cost  of  the  accelerator  and  the 
housing  facility  may  result  if  high-current,  compact  accel¬ 
erators  could  be  used  to  generate  the  high-energy  proton 
beam.  The  potential  of  high-current  accelerator  technology 
to  this  important  application  is  presently  under  evaluation. 
Results  will  be  reported  in  a  future  publication. 

In  conclusion,  accelerators  have  several  interesting  in¬ 
dustrial  applications.  However,  the  majority  of  these  appli¬ 
cations  require  low  unit  irradiation  cost.  High-current,  com¬ 
pact  accelerators  have  the  potential  to  reduce  the  unit 
irradiation  cost,  and  therefore  their  development  is  a  step  in 
the  right  direction. 
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The  experimenuil  results  on  the  trapping  of  the  beam  in  the  Naval  Research  Laboratory  modified  be¬ 
tatron  accelerator  are  in  good  agreement  with  a  revised  model  of  resistive  trapping,  and  thus  it  may  be 
concluded  that  the _^wall  resistivity  is  responsible  for  the  inward  spiral  motion  of  the  beam  after  injection. 


INTRODUCTION 

Currently,  several  laboratories  [1-5]  are  engaged  in 
studies  to  assess  the  feasibility  of  compact,  high-current 
accelerators.  Among  the  various  accelerating  schemes 
presently  under  investigation  is  the  modified  betatron  ac¬ 
celerator  (MBA)  [6-8].  This  device  is  under  study  at  the 
University  of  California,  Irvine  [2]  and  also  at  the  Naval 
Research  Laboratory  (NRL).  Since  the  initial  successful 
demonstration  of  acceleration  [1]  approximately  two 
years  ago,  the  NRL-MBA  has  achieved  [11]  trapped 
currents  as  high  as  1.5  k A  and  energies  approximately  18 
MeV.  The  beam  lifetime  that  is  typically  100-900  fisec  is 
limited  by  the  cyclotron  resonance. 

Following  the  installation  of  strong  focusing  windings 
[9,10]  in  the  NRL  device  it  is  routinely  observed  [1,11] 
that  for  several  combinations  of  injection  parameters  the 
beam  consistently  spirals  from  the  injection  position  to 
the  magnetic  minor  axis  and  is  trapped.  The  explanation 
of  this  interesting  phenomenon  has  been  so  far  elusive. 
However,  a  fair  understanding  of  the  trapping  mecha¬ 
nism  is  not  only  of  academic  interest  but  a  necessity  for 
any  upgrading  of  the  existing  or  the  construction  of  a 
new  device. 

In  this  paper  we  report  recent  experimental  results  on 
the  trapping  of  the  beam  in  the  NRL-MBA.  The  results 
are  in  agreement  with  a  revised  model  of  resistive  trap¬ 
ping  [12].  Two  modifications  have  been  introducH  to  the 
original  model.  First,  the  beam  motion  is  not  limited 
near  the  minor  axis  and  therefore  nonlinear  effects  and 
the  fast  diffusion  times  that  scale  as  —a)^ /v^p, 

where  b  —  a  is  the  thickness  of  the  chamber  and  p  is  the 
wall  resistivity,  become  important.  Second,  in  order  to 
take  into  account  the  intermediate  motion  [10]  of  the 
beam  that  has  been  omitted  in  the  calculation  of  the  im¬ 
age  fields  of  the  beam,  the  wall  surface  resistivity  is  com¬ 
puted  using  the  skin  depth  that  corresponds  to  the  fre¬ 
quency  of  the  intermediate  mode  and  not  the  actual 
thickness  of  the  chamber. 

There  are  three  distinct  groups  of  diffusion  times  with 
which  the  self-magnetic -field  of  the  beam  leaks  out  of  a 
resistive  torus.  The  shortest  are  the  “plane”  characteris¬ 
tic  times 

'r„k -a)  /tt  pk  -  ^  , 


where  /c  =  1,2,3, ....  The  “plane”  diffusion  times  are 
important  when  the  beam  is  near  Ihe  wall.  The 
“cylinder”  diffusion  times  are 

=T^ /m=p(fi{b  ~a)/2pm  , 

where  m  =  l,2,3,  ...  .  Both  and  determine  the 
speed  with  which  the  self-magnetic-field  of  the  beam 
penetrates  the  wall  of  the  chamber  and  are  instrumental 
in  the  resistive  trapping  of  the  beam.  Finally,  the  “loop” 
diffusion  time  Too=t,  =  2Tj.(ln(8ro/a )  — 2],  where  Tq  is 
the  major  radius  of  the  torus,  determines  the  speed  with 
which  the  beam  field  diffuses  into  the  hole  of  the 
doughnut.  The  “loop”  diffusion  time  docs  not  play  any 
role  in  the  resistive  trapping  of  the  beam. 

The  resistive  trapping  is  due  to  the  negative  radial 
component  of  the  image  magnetic  field  of  the  beam  that 
acts  on  its  centroid,  when  such  a  beam  moves  poloidally 
inside  a  resistive  chamber.  This  field  component  crossed 
with  the  axial  (toroidal)  velocity  of  the  beam  produces  a 
poloidal  force,  which  is  in  the  opposite  direction  to  the 
poloidal  motion  of  the  beam.  In  the  absence  of  the  strong 
focusing  and  when  the  self-fields  dominate  the  external 
fields  (high-current  regime),  the  poloidal  force  in  conjunc¬ 
tion  with  the  axial  (toroidal)  magnetic  field  drives  the 
beam  to  the  wall  (drag  instability  [13]).  However,  in  the 
presence  of  strong  focusing  the  direction  of  the  poloidal 
motion  can  be  reversed  and  Ihe  beam  spirals  to  the  minor 
axis  [12]. 

BRIEF  DESCRIPTION  OF  THE  EXPERIMENT 

The  NRL  modified  betatron  has  been  described  [1,14] 
previously.  In  this  paper  we  give,  for  completeness,  a 
short  description  of  its  basic  components.  The  NRL- 
MBA  is  a  toroidal  device  that  comprises  three  different 
external  magnetic  fields;  the  betatron  field  B,  that  can 
vary  from  0-2.7  kG,  the  toroidal  field  that  can  vary 
between  0-5. 1  kG,  and  the  strong  focusing  field  that  has 
a  maximum  gradient  between  0-31  G/cm,  when  the 
current  /jp  in  the  windings  varies  from  0-30  kA. 

The  l(X)-cm  major  radius,  I5.2-cm-inside  minor  radius 
vacuum  chamber  has  been  constructed  using  epoxy- 
reinforced  carbon  fibers.  The  desired  conductivity  is  ob¬ 
tained  by  embedding  in  the  outer  layer  of  graphite  a 
phosphor  bronze  screen.  The  measured  dc  resistance  of 
the  toroidal  vessel  is  68  ±2  mfl  and  the  corresponding 
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surface  resistivity  is  10.3  mfl  on  a  square.  The  graphite 
surface  resistivity  is  26.6  mfl  on  a  square. 

The  electrons  are  emitted  from  one  end-face  of  a  cylin¬ 
drical  carbon  cathode.  The  other  end-face  is  mounded  on 
the  cathode  stalk.  The  emitting  surface  of  the  cathode 
faces  the  circular  opening  of  the  conical  anode,  which  is 
located  on  the  midplane  of  the  device,  8.7  cm  from  its 
minor  axis. 

EXPERIMENTAL  RESULTS 

Over  a  wide  range  of  parameters  and  after  fine  tuning 
the  external  fields  the  beam  -pirals  from  the  injector  near 
the  minor  axis  and  is  trapped.  The  beam  trapping  time, 
i.e.,  the  time  it  takes  the  beam  to  travel  from  the  injection 
position  to  the  vicinity  of  the  minor  axis  is  determined  by 
measuring  the  time  delay  between  the  x-ray  peaks  that 
are  generated  at  injection  and  at  a  IXl.l-cm,  0.8-mm- 
thick  lead  target  that  is  located  on  the  magnetic  minor 
axis.  The  lead  target  is  mounted  on  the  front  surface  of  a 
3-/Ltm-thick  polycarbonate  foil  that  is  stretched  across  the 
minor  cross  section  of  the  vacuum  chamber  as  snown  in 
Fig.  1(a).  The  x  rays  arc  monitored  by  a  collimated  x-ray 
detector  that  is  located  4  m  away  from  the  lead  target. 
The  scintillator-photomultiplier  tube  is  housed  inside  a 
lead  box  and  the  x  rays  enter  the  scintillator  through  a  3- 
mm-diam  hole.  The  foil  is  graphite  coated  on  the 
upstream  side  to  avoid  charging.  Figure  Kb)  shows  an 
open-shutter  photograph  of  the  light  emitted  as  the  beam 
passes  through  the  foil.  The  x  rays  emitted  as  the  beam 


strikes  the  diode  and  the  lead  target  are  shown  in  the 
upper  trace  of  Fig.  1(c).  The  trapping  time  for  this 
shot  is  1.25  fisec.  The  lower  trace  of  Fig.  1(c)  shows  the 
output  of  the  Rogowski  coil  that  monitors  the  beam 
current.  The  peak  of  the  signal  corresponds  to  1.2  IcA. 

The  results  shown  in  Fig.  1  were  taken  with  a  0.5-cm 
hole  in  the  anode.  This  hole  is  by  a  factor  of  3  smaller 
than  that  used  regularly  in  the  NRL  device.  Thus,  the 
trapped  current  has  been  reduced  by  a  factor  of  2-3. 
This  reduction  in  the  beam  current  was  necessary  in  or¬ 
der  to  achieve  satisfactory  resolution  in  the  open-shutter 
photographs. 

To  determine  the  effect  of  the  foil  on  the  transverse 
beam  orbit,  we  carried  out  a  series  of  experiments  in 
which  the  3-^im-thick  foil  was  replaced  with  a  foil  of  the 
same  composition  but  with  only  half  its  thickness.  The 
results  show  that  the  equilibrium  position  of  the  beam  is 
slightly  larger  in  the  case  of  1.5-/im-thick  foil.  It  requires 
approximately  1-2  G  higher  vertical  field  (~4-8%)  to 
shift  the  equilibrium  to  its  original  position  and  make  the 
orbits  identical. 

As  the  electrons  pass  through  the  plastic  foil,  they 
suffer  both  inelastic  and  elastic  scattering.  The  stopping 
power  [15]  of  0.6  MeV  electrons  passing  through  po¬ 
lyethylene  (data  for  polycarbonate  are  not  available)  is 
~2  (MeV  cm^)/gm.  Thus,  the  energy  loss  per  pass  is  0.6 
keV.  The  total  energy  loss  in  1.2  /isec,  i.e.,  in  sixty  revo¬ 
lutions  around  the  major  axis,  is  36  keV  or  ~6%.  The 
energy  loss  in  the  thinner  foil  is  only  18  keV  and  there¬ 
fore  the  equilibrium  position  is  expected  to  increase  by 
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BtiVr.) 

FIG.  2.  (a)  Trapping  time  reproducibility  and  (b)  dependence 
of  T„  and  Tj  on  Bg  for  constant  /Bg. 


3%  when  the  thickness  of  the  foil  is  reduced  to  half.  This 
shift  is  not  substantially  different  from  that  observed  in 
the  experiment.  The  elastic>scattering  induced  RMS  an¬ 
gle  is  0.9*  for  the  first  pass  through  the  KS-pm-thick  foil. 
Although  substantial,  the  elastic  scattering  does  not  con¬ 
tribute  to  the  shift  of  the  equilibrium  position. 

The  beam  orbits  are  very  reproducible  and  shows 
only  modest  variations  for  the  same  operating  parame¬ 
ters.  Figure  2(a)  shows  for  seven  shots  taken  with  the 
same  values  of  the  fields.  It  is  apparent  from  this  figure 
that  T„  varies  by  ±7%.  In  a  second  run  with  five  shots 
the  variation  was  even  smaller.  Figure  2(bi  shows  t„,  vs 
Bg  for  constant  I^f/Bg,  where  is  the  current  of  the 
strong-focusing  windings.  For  all  practical  purposes 
remains  constant  in  the  narrow  range  tested. 

In  addition  to  t„,  the  bounce  period  Tg,  i.e.,  the  time 
the  beam  takes  to  perform  a  complete  revolution  in  the 
poloidal  direction,  is  of  special  interest.  To  determine  Tg, 
the  foil  target  was  removed  and  replaced  with  a  1.1-cm- 
wide,  1-mm-thick,  16-cm-long  lead  strip.  The  lead  target 
is  backed  on  the  upstream  side  by  a  thin  plastic  strip  and 
is  mounted  on  a  half  lucite  ring  that  is  carbon  coated. 
The  symmetry  axis  of  the  target  lies  on  the  midplane  of 
the  device  as  shown  in  Fig.  3(a). 

The  light  emitted  from  the  upstream  side  of  the  target 
when  the  beam  strikes  it  is  monitored  with  an  open 
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shutter  camera.  Results  arc  shown  in  Fig.  3(b).  The  x- 
ray  signal  and  the  output  of  the  Rogowski  coil  that  moni¬ 
tors  the  beam  current  are  shown  in  Fig.  3(c).  The  bounce 
period  is  inferred  from  the  time  delay  of  the  two  x-ray 
peaks,  as  indicated  in  Fig.  3(c),  and  in  this  shot  is  840 
nsec.  The  damage  pattern  on  the  lead  strip  has  a  diameter 
that  is  equal  to  the  diameter  of  the  anode  hole.  This  im¬ 
plies  that  the  diameter  of  the  beam  has  not  changed  after 
about  40  revolutions  around  the  major  axis.  In  addition, 
we  observe  that  the  damage  pattern  is  a  semicircle  that  is 
always  located  near  the  lower  edge  of  the  strip.  From 
this  observation  it  may  be  concluded  that  the  beam  drifts 
3  mm  over  20  nsec,  i.e.,  its  bounce  speed  near  the  strip  is 
~  15  cm//isec. 

To  verify  that  there  is  no  correlation  between  r,,  and 
Tf,  i.e.,  with  the  speed  the  beam  magnetic  field  diffuses 
into  the  hole  of  the  doughnut,  the  vacuum  chamber  was 
unboiled  in  two  joints  that  are  located  180“  apart  in  the 
toroidal  direction  and  a  ring  insulator  was  inserted  in 
each  of  these  joints.  Sixty  carbon  resistors,  51  ft  each, 
were  symmetrically  mounted  on  the  outer  surface  of  one 
of  the  two  rings  as  shown  in  Fig.  4(a).  To  improve  its 
voltage  holding  capabilities  the  inner  surface  of  the  blue 
nylon  insulator  was  angled  and  a  0.6-cm-deep  groove  was 
machined  at  its  plane  of  symmetry.  In  addition,  its  inner 
surface  was  protected  from  stray  electrons  by  a  0.8-mm- 
thick  lead  strap  that  is  supported  by  an  epoxy  reinforced 


carbon  fiber  belt.  The  purpose  of  the  second  insulator 
was  only  to  minimize  the  distortion  of  the  toroidal 
chamber  and  thus  shorting  wide  straps,  instead  of  resis¬ 
tors,  were  installed  on  its  outer  surface. 

There  are  two  distinct  currents  flowing  on  the  wall  of 
the  vacuum  chamber.  The  first  is  due  to  the  rising 
vertical  field  and  the  second  I,  is  induced  by  the  beam. 
Since  the  vertical  field  during  the  first  quarter  period 
varies  as  B^(f)  =  fi^sin(2rrt /t),  where  is  the  peak 
field  and  t  is  the  period,  the  induced  voltage  in  the 
chamber  is  V= -V^cosa-nt /t),  where 
t)B^,  ^nd  Tq  ts  major  radius  of  the  torus.  The 
current  flowing  on  the  wall  of  the  chamber  is  described 
by  the  equation  V  =  L{di^/dt  )  +  where  L  is  the  in¬ 
ductance  and  R  the  resistance  of  the  torus.  The  instan¬ 
taneous  value  of  can  be  found  by  integrating  the  above 
equaton  and  is  given  by 


R 

1  + 

2ttL  /R 

2 

r 

cos(2irf  /t)+  sin(2xrr/r)  —  ^ 


(ft 
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FIG  4  'ai  Schematic  of  the  vacuum  chamber  joint  with  the  insulator  in  place  fb'  Current  and  gap  voltage  waveforms 
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with  the  initial  condition  =0)  =  0. 

The  temporal  profile  of  predicted  by  the  above  sim¬ 
ple  model  is  identical  to  that  predicted  from  the  exact 
solution  of  the  diffusion  problem  for  a  toroidal  resistive 
shell  [16]  and  also  is  in  good  agreement  with  the  results 
of  the  TRIDIF  code  for  a  finite  thickness  toroidal  vessel. 
According  to  Eq.  (1),  i„,  has  a  maximum  at  time  which 
for  a  —  2-rriL /R)/t  « I  is  determined  from 
a^[  1 —r^ /(£.//?)]= —e  of  the 

current  is  i^=y/R.  The  measured  peak  value  of  the 
current  in  the  experiment  is  in  good  agreement  with  the 
above  theoretical  prediction  and  scales,  as  expected,  with 
the  value  of  the  resistance  at  the  gap. 

The  return  current  of  the  beam  is  measured  with  a  fast 
Rogowski  coil  (7,^5^  =  20  nsec)  that  is  located  outside  the 
vacuum  chamber.  With  the  resistors  at  the  gap  shorted, 
the  Rogowski  coil  shows  a  slowly  rising  current  that  is 
consistent  with  the  decay  of  /,.  However,  when  the 
shorting  clips  are  removed  the  Rogowski  coil  shows  a 
current  pulse  that  rises  to  of  its  peak  value  in  less 
than  100  nsec  as  shown  in  Fig.  4(b).  The  lower  trace  in 
fig.  4(b)  shows  the  voltage  across  the  resistors  Vg  as  mea¬ 
sured  directly  by  a  Tektronix  7844  oscilloscope  after  a 
1(X)A’  attenuation.  The  shape  of  the  time-integrated  Vg  is 
very  similar  to  the  current  waveform  resistered  by  the 
Rogowski  coil,  i.e.,  Vg  is  proportional  to  the  time  deriva¬ 
tive  of  the  current. 


COMPARISO.N  WITH  THEORY' 

The  equations  that  describe  the  motion  of  the  beam 
centroid  base  been  solved  numerically  using  analytical 
expressions  for  R.  and  fig.  The  stellarator  fields  are 
determined  numerically  from  Biot-Savart  law  by  dividing 
each  period  of  the  windings  into  20  segments.  The  image 
fields  on  the  beam  centroid  have  been  computed  analyti¬ 
cally  for  a  uniform  density  electron  ring  that  is  located 
inside  a  large  major  radius  torus  with  resistive  wall  of 
thickness  b-a,  where  a  is  the  inner  and  b  the  outer  radii  of 
the  torus.  In  contrast  with  previous  calculation  [12],  the 
beam  is  not  limited  near  the  minor  axis. 

In  the  local  cylindrical  coordinate  system  (p.6.z)  with 
its  origin  on  the  geometric  minor  axis,  the  electrostatic 
potential  inside  the  ring  is  given  'ly  [17] 
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where  A  and  a  define  the  beam  position  on  the  transverse 
plane,  is  the  minor  radius  of  the  beam,  and  Qi  is  the 
charge  per  unit  length.  Similarly,  the  magnetic  vector 
potential  inside  the  beam  is  [17] 
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where  Rg  is  the  normalized  toroidal  beam  velocity. 

The  time-dependent  coefficients  are 

zero  at  r  =0  and  are  determined  by  the  differential  equa¬ 
tions 
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where 
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a  is  the  wall  conductivity  of  the  toroidal  chamber,  and 
are  the  zeros  of  the  function 
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when  m  =0.  and 
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when  ,t7  =  1,2.3 .  The  arguments  of  Bessel  func¬ 

tions  in  (5)  and  (61  are  x^=aa,  x,  —ab  and  fip  is  the  ma¬ 
jor  radius  of  the  beam.  For  each  m  there  is  an  infinite 
number  of  zeros  denoted  by  the  index  A;  =0,1.2  The 
time-dependent  coefficients  that  appear  in  the 

vector  potential  are  equal  to 
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where /,[,(a)  is  the  dematise  of  f^^a)  and 
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when  m  =  1,2, 3 . The  image  fields  at  the  beam  cen¬ 

troid,  i.e.,  when  p=A  and  ^=a,  are  obtained  from  the 
expressions  of  4>o  and  A  ^  given  above. 

The  radial  components  at  the  beam  centroid  is  of 
special  interest  because  it  is  responsible  for  the  inward  ra¬ 
dial  motion  of  the  beam.  This  component  is  given  by 


i  i  ^ 

”  a  a 


X  [  —  Ul^l{t)cosma{t) 

+  Ul^lithinmait)]  . 


TABLE  1.  Parameters  of  the  run  shown  in  Fig,  i 

Toms  major  radius  100  cm 

Torus  minor  radius  a  15.2  cm 

ReUdvistic  factor  y  1.5 

SF  radius  pg  23.4  cm 

SF  current  J^f  24  kA 

Vertical  field  at  injection  26  G 

Toroidal  field  Bgo  4  kG 

Beam  minor  radius  3  mm 

Beam  current  /»  1.2  kA 

Wall  resistivity  8  mflcm 

Intermediate  frequency,  _ 1.8X  10’  sec' 


Since  ’  is  independent  of  and  C/o*  ,  i.e.,  it  is  in-  To  gain  further  insight,  we  have  computed  the  image 
dependent  of  loop  time  t,  =  Too,  the  beam  trapping  time  fields,  including  first-order  toroidal  corrections,  at  the 
should  also  be  independent  of  t,  as  observed  in  the  exper-  beam  centroid  for  a  beam  inside  a  resistive  toroidal 
iment.  shell — in  this  case 
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where  the  time-dependent  parameters  U'qHi),  f/g'lr), 
O'),  O*').  O'),  and  0(/)  are 

zero  at  /  =  0  and  are  determined  by  the  differential  equa¬ 
tions 


1  4_ 

1  .1.  ‘ 

°  AttcqC  t, 

1  ^ 

2Rg 

+  _L  //'  1 '  = _ !_  f/io< 


,-,ioi  ,  rn  _  ^QiPo  m  A 

+  —  f-'en.  -  - — 

^v(qC  a 

lj,0'  I  w  rn  A 


47760^  T, 


ff'  I  ’4.  f'  !  —  _ 

^rm  o  ^cm 

i.  IKh 


O'"  + 

'’’f 

'  I '  +  Z!L  V  I '  =  -±- 

Ifin-foC  r,  R^, 


cost  m  —  1  )a 


1/ 1 1 1  .  i/ii >  = 

^  jm  ~  ^  ^  jm 


QiPd  m  ^  _a_  ^ 
16-foC  Rt,  a 

[  m  —  1 

A 

X  —  sin(m  —  1  )a  , 

o 


where 
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x^=Acosa  =  Bj,  — Tq  . 

It  is  apparent  from  Eqs.  (8)  and  (9)  that  the  toroidal 
correction  term  of  is  a  function  of  t,.  However,  this 
term  is  multiplied  by  sina(r)  and  therefore  averages  to 
zero  in  a  poloidal  period. 

In  the  limit  (b  —  fl)«fl,  the  toroidal  electric  field  £p 
on  the  outer  surface  of  the  chamber  for  a  stationary  beam 
has  a  relatively  simple,  closed  form.  At  f  =0,  £@  =  0  and 
peaks  within  a  few  fast  diffusion  times  t^.  For  a  longer 
time.  Eg  decays  to  zero  with  t,.  This  form  of  the  electric 
field  is  consistent  with  the  observed  return  current  after 
the  beam  injection.  When  the  resistors  at  the  gap  are 
shorted  the  current  measured  with  the  external  Rogowski 
coil  should  rise  at  the  same  rate  the  beam  return  current 
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FIG.  5.  Beam  centroid  orbit  from  the  numerical  integration 
of  the  equations  of  motion,  using  the  image  fields  from  the  resis¬ 
tive  shell  model  [(a)  and  (b)].  Results  from  the  experiment  (c). 

decays.  However,  when  the  resistors  at  the  gap  are  not 
shorted,  a  portion  i/,,)  of  the  return  current  (/,)  flows  on 
the  outer  wall  of  the  chamber  as  shown  in  Fig.  4.  rises 
considerably  faster  because  the  beam  magnetic  field  can 
leak  locally  out  of  the  resistive  gap  considerably  faster 
than  it  can  leak  out  of  a  uniform  chamber. 

Results  from  the  numerical  integration  of  orbit  equa¬ 
tions  are  shown  in  Fig.  5.  The  various  parameters  for  the 
run  are  listed  in  Table  I.  Figure  5(a)  shows  the  projection 
of  the  centroid’s  orbit  on  the  6  =  0  plans.  Both  the  inter¬ 
mediate  and  slow  (bounce)  modes  are  apparent.  Since 
there  are  six  field  periods  between  O’SO’Slir,  the  elec¬ 
trons  perform  six  oscillations  during  one  revolution 
around  the  major  axis.  To  take  into  account  the  inter¬ 
mediate  motion  that  has  been  neglected  in  the  calculatoin 
of  the  image  fields,  the  surface  resistivity  in  the  code  is 


computed  using  the  skin  depth  that  corresponds  to  the 
intermediate  frequency  and  not  the  actual  thickness  of 
the  wall. 

The  solid  circles  in  Fig.  5(b)  show  the  positions  the 
beam  crosses  the  0  =  240*  plane.  This  is  a  realistic  simu¬ 
lation  of  the  experimental  situation.  The  time  difference 
between  two  circles  is  equal  to  the  period  around  the  ma¬ 
jor  axis,  i.e.,  ~23  nsec.  The  parameters  of  this  run  are 
similar  to  those  in  Fig.  5(c)  and  the  similarity  of  the  two 
orbits  is  quite  apparent.  When  the  crossing  plane  is 
moved  from  0=240'  to  a  different  azimuthal  position  0, 
the  beam  orbit  rotates  around  the  minor  axis.  The  rota¬ 
tion  predicted  by  the  theory  is  very  similar  to  that  ob¬ 
served  in  the  experiment. 

There  is  some  ambiguity,  both  in  the  experiment  and 
theory,  in  the  determination  of  the  beam  trapping  time, 
because  its  exact  value  depends  on  the  position  and  size 
of  the  target.  However,  this  is  not  the  case  with  xg, 
which  can  be  measured  very  accurately.  We  made  four 
computer  runs  for  different  values  of  Bg  keeping 
/sF/iJe  =  constant.  Figure  2(b)  shows  Xg  vs  Bg  for  three 
of  these  runs.  For  all  practical  purposes  Xg  remains  con¬ 
stant  as  Rg  varies.  In  the  fourth  run  Bg  was  reduced  to  2 
kG  and  although  the  beam  orbit  changed  substantially  Xg 
was  lower  only  by  7%. 

As  a  rule,  the  theory  predicts  a  Xg  and  that  are  ap¬ 
proximately  a  factor  of  2  shorter  than  those  observed  in 
the  experiment.  With  the  exception  of  these  two  times 
the  revised  model  of  resistive  trapping  is  in  agreement 
with  the  experiment  observations,  although  in  the 
analysis  the  beam  current  remains  constant  while  in  the 
experiment  the  current  decays.  This  decay  is  observed 
whenever  there  is  a  target  inside  the  chamber  but  in  gen¬ 
eral  is  absent  during  the  acceleration  experiments  when 
the  various  targets  are  removed. 
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Excitation  of  the  t  —  VI  Cyclotron  Resonance  in  the  NRL 
Modified  Betatron  Accelerator* 

L.  K.  Len,*^  T.  Smith, P.  Loschialpo,  J.  Mathew,  S.  J.  Marsh,^^  D.  Dialetis,') 

J.  Golden,*^ ^  J.H.  Chang  and  C.  A.  Kapetanakos 

Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  DC  20375-5000 

The  basic  periodicity  in  the  NRL  device  is  twelve-fold.  There  are  twelve  toroidal 
field  coils,  twelve  sectors,  and  so  on.  For  a  long  time,  the  £  =  12  resonance 
was  the  dominaint  resonance  in  the  experiment.  To  test  the  importance  of  the 
£  =  12  resonance,  we  have  intentionally  introduced  an  £  =  12  field  error  using 
twelve  resonant  coils.  By  activating  these  coils  when  the  resonance  is  observed, 
the  duration  of  x-rays  produced  by  beam  loss  was  reduced  from  900  /is  to  5.5 
ps,  while  the  amplitude  of  the  signal  increased  from  0.5  to  ~40  volts.  The  full 
width  of  the  x-ray  pulse  at  half  maximum  is  inversely  proportional  to  the  current 
through  the  resonant  coils  and  proportional  to  the  risetime  of  the  pulse.  Work  is 
in  progress  with  a  set  of  twelve  internal  coils  that  have  a  risetime  of  300  -  400 
ns.  The  results  with  these  coils  have  confirmed  the  importance  of  the  £  =  12 
resonance.  In  addition,  it  has  been  shown  experimenteJly  that  by  adding  a  ABg 
such  that  the  ratio  of  the  total  toroidal  magnetic  field  to  the  vertical  field  is 
not  an  integer,  the  cyclotron  resonance  is  not  excited.  These  experiments  have 
extended  the  beam  lifetime  by  more  than  100  fxs,  which  is  approximately  equal 
to  the  risetime  of  the  applied  ABg  pulse.  In  the  presence  of  such  a  pulse,  the 
beam  lifetime  is  1  ms  and  the  electron  beam  energy  is  22  MeW  while  the  trapped 
current  is  in  excess  of  1  kA. 

I.  INTRODUCTION 

The  lifetime  of  the  electron  beam  in  the  NRL  modified  betatron  accelerator  has  been 
limited  by  the  excitation  of  the  cyclotron  resonances.’  These  resonances  occur  whenever 
the  condition  (BgoJ B:o )  =  (2£^  -  l)/2£  is  satisfied,  with  (  taking  integer  values. Here  Bgo 
and  Bzo  are  the  toroidal  and  vortical  magnetic  fields  on  the  minor  axis  of  the  chamber.  In 
this  paper,  we  present  results  from  several  experiments  related  to  the  excitation  of  electron 
cyclotron  resonances.  The  first  experiment  deals  specifically  with  the  £  =  12  resonance  and 
the  .second  addresses  the  spatial  distribution  of  the  beam  los.ses.  In  the  third  experiment, 
a  time  rising  ABg  has  been  added  to  the  main  Bg  field.  This  has  effectively  resulted  in  a 
total  Bgr  that  maintains  the  condition  of  Bg-f/B^  =  constant  ^  integer  over  the  risetime 
of  the  ABg  pulse,  and  thus  the  beam  does  not  cross  any  resonance. 

II.  EXPERIMENTAL  APPARATUS 

Detailed  descri]>tion  of  the  NRL  modified  betatron  accelerator  can  be  found  in  pre¬ 
vious  publications’ '  where  various  cxi)erimental  results,  incl\iding  ol)servation  of  the 
electron  cyclotron  resonances  were  rejiorted.  There  are  thn'e  main  external  magnetic 
fields  that  jjrovide  confinement  to  the  high  current  beam,  namely,  the  toroidal  (Bg).  the 
strong  focusing  (£?.«/).  and  the  vertical  (B-)  field.s.  They  are  generated  by  pulsed,  aircore 
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electromagnets  and  have  risetimes  equal  to  2.3  ms,  1.0  ms,  and  2.6  ms  respectively.  The 
peak  value  of  Be  is  ~5  kG  and  that  of  B,  is  ~2  kG.  The  strong  focusing  windings  can  be 
operated  at  peeik  currents  ais  high  as  ~30  kA. 

Two  sets  of  resonant  coils  have  been  used  in  the  experiments  reported  in  this  paper. 
A  set  of  12  single-turn  coils  are  installed  on  the  outside  of  the  vacuum  chamber  (external 
resonant  coils).  They  are  approximately  equally  spaced  around  the  toroidal  chamber.  The 
drivers  for  these  coils  can  be  configured  to  generate  ABe  with  risetime  of  12  and  5  ps.  To 
generate  a  faster  resonant  pulse,  these  coils  have  been  replaced  with  12  single-turn  internal 
coils  that  are  located  accurately  at  the  joints  between  sectors  as  shown  in  Fig.  1.  The 
internal  coils  are  wound  on  a  blue  nylon  housing  and  encapsulated  with  epoxy. 


FRj.  1  Sch<n\<  ti<  \  it-u  ol  rrj^onatit  roils  with  rosp<‘rt  !o  the  vacuum  citHinhrr  spc^jors 


In  tli('  <'X])erinienf  with  B^rlB-  —  constant  integer,  the  nunilier  of  external  coils 
wa.s  increas('(l  to  24  in  order  to  jiroduce  a  smoother  ABe  In  addition,  the  pulsed  power 
driver  was  modified  to  inrrea.se  the  quarter  period  risetiine  to  "^lOO  /iS. 
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Beam  injection  is  initiated  near  the  peak  of  the  toroidal  and  strong  focusing  fields. 
The  injector  diode  is  located  inside  the  vacuum  chamber  8.6  cm  from  the  minor  axis. 
When  the  ratio  of  Bg  to  is  nearly  12,  the  resonant  coils  are  energized  to  produce  the 
ABg  disturbance.  Diagnostics  includes  the  following:  Magnetic  field  probes  to  measure 
Be  and  B~,  three  NPM-54  x-ray  detectors  that  monitor  the  x-rays  generated  when  the 
beam  strikes  the  chamber  wall,  a  Rogowski  belt  to  measure  the  electron  beam  current  and 
optical  fibers  that  are  located  on  the  outside  of  the  vacuum  chamber.  The  light  emitted 
when  the  energetic  electrons  strike  these  fibers  is  monitored  with  photomultiplier  tubes. 
This  technique  has  allowed  us  to  determine  accurately  the  beam  dump  positions. 

III.  EXPERIMENTAL  RESULTS 

In  the  modified  betatron  the  electron  beam  performs  a  complicated  motion  following 
injection.  To  avoid  hitting  the  diode  or/and  the  wall  as  the  beam  spirals  onto  the  minor 
axis  during  the  trapjring  process,  it  is  necessary  to  select  the  operating  parameters  carefully 
to  maximize  the  beam  trapping  efficiency.  The  signals  from  the  x-ray  detectors  looking  at 
the  diode  and  the  lead  shielding  around  the  Rogowski  belt  are  used  in  conjunction  with 
the  beam  current  monitor  as  a  guide  to  arrive  at  the  optimum  field  settings.  The  beam 
lifetime  is  inferred  from  the  x-ray  signal  such  as  that  shown  in  Fig.  2.  The  main  peak  at 
800  //s  corresironds  to  an  energ^v  of  20  Me\'  while  the  last  peak  at  1300  //s  is  28  Me\h 


Tim<>  (fin) 

FK  'f  2  I  v|)iral  x-ray  signal 


For  a  long  time,  the  (  --  12  resonance  has  been  the  dominatit  cause  of  beam  loss.  Even 
after  a  careful  realignment  of  the  fiehl  coils,  vacuum  chamber,  the  injector,  the  sui)porting 
deck  strtictures  and  nahiction"-  of  other  possible  field  errors,  tlu'  (  ~  12  n'sonance  remain- 
significant.  It  i-  t  h<'  fiist  major  resonam'e  ('ncounf<'red  as  can  be  seen  in  Fig.  2. 

The  (  =  12  resonauc('  is  of  specitd  int<-rest  because  the  XRL  cU'vice  has  a  twelve-fohl 
I)eriofhcit y.  There  ate  twehe  s<'cfors  in  tin-  vacuum  chamber  ;uid  tw(’h('  toroifhd  field  coils. 
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Time  (fis) 

FIG.  3,  X-raj  .'iignal  wlipn  thp  resonani  coils  are  activated  at  <'=12.  (a)  Full  trace,  (b)  Expanded 
trace. 


To  enhance  this  resonance,  wc  have  installed  twelve  resonant  coils  to  generate  a  periodic 
field  perturbation. 

Initially,  a  set  of  twelve  resonant  coils  was  installed  on  the  outside  of  the  chamber 
(external  coils).  These  coils  are  initiated  at  r  =  430  //s.  i.e,.  when  the  f  =  12  condition  is 
reached.  A  ADe  pulse  with  a  quarter  period  risetiine  of  12  ps  and  amplitude  of  200  G  is 
generated  on  the  minor  axis  when  the  current  in  the  coils  is  9  kA.  As  it  is  apparent  from 
Figs.  2  and  3(a).  the  width  of  the  .x-ray  signal  (a  mea.sure  of  rate  of  beam  loss)  is  reduced 
from  approximately  900  //s  to  S  /is.  i.e..  by  more  than  two  orders  of  magnitude  while  its 
amplitude  has  increase  by  a  factor  of  twenty.  The  width  of  the  x-ray  signal  varies  inversely 
to  the  current  in  the  resonmit  coils  as  .shown  in  Fig.  4.  The  loss  rate  is  also  dependent  on 
the  risetime  of  the  resonant  coil  pulsix 

By  connecting  the  coils  in  parallel  instead  of  in  series,  the  risetime  was  reduced  to  5 
/IS.  and  even  shorter  risetime  pulses  have  been  produced  with  a  set  of  internal  coils.  The 
internal  coils  are  powered  by  new  drivers  with  a  risetime  of  400  ns.  Being  faster  and  nearer 
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Resonant  Coil  Current  (kA) 

FIG.  4.  Full  width  of  the  x-ray  pulse  at  half-maximum  versus  the  current  in  the  resonant  coils.  The 
peak  ^3$  generated  in  the  chamber  corresponds  to  28  G/kA. 


RC  Current  Risetime  (/is) 

FIG.  5.  Full  widili  of  tlic  x-ray  pulse  at  half  maximum  versus  the  current  riselime  of  the  resonant 
coils. 


to  the  sector  flanges,  the  internal  coils  require  substantially  higher  voltage  and  current  to 
produce  tlic  required  field  perturbations.  The  results  from  this  experiment  are  shown  in 
Fig.  5. 

To  determine  the  spatial  distribution  of  beam  losses  when  the  resonant  coils  are  ener¬ 
gized.  several  400-//ni-dlameter  optical  fibers  were  mounted  on  the  outside  of  the  vacuum 
chamber.  By  the  time  the  f  =  12  resonance  is  crossed,  the  electrons  have  acquired  suf¬ 
ficient  energy  to  penetrate  the  chamber.  The  light  generated  when  the  electrons  strike 
the  fiber  is  inouitoif'd  with  a  photomultiplier  tube.  The  results  from  scanning  around  the 
torus  at  10°  intervals  are  shown  in  Fig.  6.  Only  six  peaks  instead  of  the  '•xpected  t.velve 
are  observed.  This  is  most  likely  due  to  the  beating  between  the  f  =  12  and  f  =  6  mode 
that  is  excited  by  tlu'  six  field  periods  of  the  stellarator  windings.  Computer  simulation 
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shows  that  the  strong  focusing  modulates  the  ^  =  12  peaks  in  such  a  way  that  only  every 
other  strikes  the  wall  as  shown  in  Fig.  7.  It  is  still  not  clear  at  this  time  why  three  of  the 
peahs  are  higher  than  the  rest.  Experiments  are  in  progress  to  clarify  this  issue. 


V»e»  Orbits 


e 

FIG.  T.  Results  from  computer  simulation  Rq  is  the  major  radius. 


Finally,  wo  pro.sont  tho  re.sults  from  an  expeiimont  in  which  a  time-rising  is  used 
to  extend  the  beam  lifetime.  It  is  well  known  that  when  Bsj  B.  —  20  l{2C—\ )  ^  C  =  integer 
(for  the  C  values  relevant  to  our  observation),  resonances  are  excited  and  beam  losses  are 
observed.  Therfore,  by  using  an  applied  ABe  to  maintain  Bot/B-  =  constant  ^  integer. 
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where  Ber  =  Bg  ABg,  crossing  of  resonances  can  be  avoided.  A  set  of  24  single-turn 
external  coils,  powered  by  a  larger  capacitor  bank  generates  the  ABg  ramp.  Results  from 
the  experiment  are  shown  in  Fig.  8  where  the  ABg  was  introduced  at  r  =  800/is.  Beam 
loss  is  suppressed  for  100  /xs,  which  is  the  risetime  of  the  applied  Bg. 


Time  (/xs) 

FIG.  8.  X-ray  signal  for  the  rase  where  a  ABg  has  been  added  to  extend  the  peak  at  830  ^ls  to  950 

liS. 
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Dynamic  behavior  of  an  electron  ring  close  to  a  cyclotron  resonance 
in  a  modified  betatron  accelerator 

D.  Dialeiis,*  S.  J.  Marsh,^  and  C.  A.  Kapctanakos 
Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  D.C.  20S7S-5000 

(Received  15  June  1992) 

The  effect  on  the  electron-ring  dynamics  when  a  cyclotron  resonance  is  crossed  in  a  modified  betatron 
accelerator  has  been  studied  analytically  and  numerically.  It  has  been  found  that,  in  the  presence  of 
small  vertical  field  errors,  there  is  a  field-error-amplitude  threshold  below  which  the  normalized  trans¬ 
verse  velocity  of  the  gyrating  electrons  is  bounded  (Fresnel  regime)  and  above  which  it  is  unbounded 
(lock-in  regime).  In  the  lock-in  regime,  the  average  value  of  the  normalized  axial  (toroidal)  momentum 
yP»,  where  y  is  the  relativistic  factor  and  pg  is  the  normalized  axial  velocity,  remains  constant,  i.c.,  the 
resonance  is  never  crossed.  In  addition,  above  threshold.  Pi  increases  proportionally  to  the  square  root 
of  the  time.  The  threshold  value  of  the  vertical  field  error  amplitude  can  be  made  larger  either  by  in¬ 
creasing  the  acceleration  rate  or  by  adding  a  small  oscillatory  toroidal  field  to  the  main  toroidal  field. 
The  multiple  crossing  of  the  same  resonance,  in  the  presence  of  such  a  small  oscillatory  toroidal  field, 
was  also  studied  with  some  interesting  results. 

PACS  number(s);  41.75.Fr,  41.85.Lc,  29.27.Bd 


I,  INTRODUCTION 

There  is  extensive  experimental  evidence  suggesting 
that  the  gradual  beam  loss  that  is  observed  in  the  Naval 
Research  Laboratory  (NRL)  modified  betatron  accelera¬ 
tor  (MBA)  is  a  consequence  of  crossing  various  cyclotron 
resonance  modes  during  acceleration  [1,2].  The  cyclo¬ 
tron  resonance  is  due  to  the  excitation  of  the  cyclotron 
motion  by  field  errors  associated  with  the  toroidal  and 
vertical  magnetic  fields.  Consequently,  these  field  errors 
can  be  either  a  vertical-field  bB.  (VF)  error  or  an  axial- 
(toroidali  field  (TF)  error  or  both. 

Recirculating  accelerators  with  low  accelerating  gra¬ 
dient  such  as  the  existing  NRL  modified  betatron  are  sen¬ 
sitive  to  field  errors,  because  the  electrons  have  to  per¬ 
form  a  large  number  of  revolutions  around  the  major  axis 
in  order  to  obtain  the  desired  peak  energy.  Successful 
detection  and  elimination  or  reduction  of  several  field  er¬ 
rors  in  the  NRL  device  led  to  beam  energies  in  excess  of 
20  MeV,  while  the  trapped  current  is  above  I  kA  [3]. 

Although  the  cyclotron  resonance  is  a  potent  mecha¬ 
nism  with  the  potential  to  disturb  the  beam  at  a  low  ac¬ 
celeration  rate  and  when  the  various  fields  are  not  care¬ 
fully  designed,  it  also  may  provide  a  powerful  technique 
for  extracting  the  beam  from  the  magnetic-field 
configuration  of  the  modified  betatron  [4],  The  study  of 
the  cyclotron  resonances  is  facilitated  by  introducing  the 
detuning  factor  w  —  rQilgQ/yPgC  —  l,  where  Cq  is  the  major 
radius  of  the  torus,  c  is  the  velocity  of  light,  y  is  the  rela¬ 
tivistic  factor.  P„  is  the  normalized  toroidal  velocity, 
n9Q=  lelBgfi/mc,  Rflf)  is  the  toroidal  magnetic  field  on  the 
minor  axis,  e  and  m  are  the  charge  and  mass  of  the  elec¬ 
tron,  and  /  is  the  mode  number  of  the  resonance.  The  / 
mode  of  the  cyclotron  resonance  occurs  when  the  ratio  of 
the  toroidal  field  5^,,  to  the  vertical  field  approxi¬ 

mately  equal  to  /.  Since  at  equilibrium  the  ratio 


raQio/yPgC  is  approximately  unity,  at  least  when  the 
beam  current  is  low  and  in  the  absence  of  strong  focus¬ 
ing,  the  /-mode  cyclotron  resonance  is  crossed  when  the 
detuning  factor  is  zero.  Furthermore,  the  detuning  factor 
appears  naturally  in  the  slow  equations  of  motion  derived 
by  averaging  out  the  fast  cyclotron  motion.  The  fact  that 
uj  +  /  is  inversely  proportional  to  yPg  has  a  profound 
effect  on  the  ring  dynamics.  The  quantity  yp^  has  a  non¬ 
linear  dependence  on  the  normalized  perpendicular  ve¬ 
locity  Pi,  and,  as  a  consequence,  there  is  a  threshold  for 
the  vertical  field  error  amplitude  below  which  P^  is 
bounded  (Fresnel  regime)  and  above  which  Pi  increases 
continuously  (lock-in  regime)  with  time.  In  the  lock-in 
regime,  the  detuning  factor  remains  almost  zero  long 
after  the  resonance  has  been  reached  and,  therefore,  the 
resonance  is  never  crossed.  In  addition,  above  threshold. 
Pi  increases  proportionally  to  the  square  root  of  time 
while  yPff  remains,  on  the  average,  a  constant.  The 
threshold  value  of  the  vertical  field  error  amplitude  can 
be  made  larger  either  by  increasing  the  acceleration  rate 
or  by  adding  a  small  oscillatory  toroidal  field  to  the  main 
toroidal  field.  The  latter  method  is  called  dynamic  stabil¬ 
ization  of  the  resonance. 

In  the  case  of  a  vertical  field  error  and  in  the  absence  of 
acceleration,  space  charge,  and  strong  focusing  field,  our 
studies  of  the  cyclotron  resonances  show  that  the  normal¬ 
ized  transverse  velocity  Pi  and  thus  the  Larmor  radius  of 
the  transverse  motion  of  the  gyrating  particles  grows 
linearly  with  time  [5j,  provided  that  nonlinear  effects  as¬ 
sociated  with  yPg  are  neglected.  When  such  effects  are 
taken  into  account.  Pi  is  periodic  and  bounded. 

As  mentioned  previously,  in  the  presence  of  an  ac¬ 
celerating  field  and  of  a  large  vertical  field  error,  P^  in¬ 
creases  proportionally  to  the  square  root  of  time  while 
yPg  saturates,  i.e.,  the  electrons  lock  into  a  specific  reso¬ 
nance  (lock-in  regime).  When  the  amplitude  of  bB.  is 


47  2043 


(c''!993  The  American  Physical  Society 


3R 


2044 


D.  DIALETIS,  S.  J.  MARSH,  AND  C.  A.  KAPETANAKOS 


47 


below  the  threshold,  0^  exhibits  Fresnel  behavior,  i.e.,  0i 
grows  quickly  for  approximately  1  /isec  and  then  satu¬ 
rates  until  the  beam  reaches  the  next  resonance. 

In  the  case  of  an  axial  field  error  and  in  the  absence  of 
acceleration,  when  the  ring  is  initially  at  exact  resonance, 
/?i  grows  exponentially  with  time  only  for  a  very  short 
period.  Since  increases  at  the  expense  of  0^  the  parti¬ 
cles  are  kicked  off  resonance  and  0^  varies  cyclicly  with 
time.  Similarly,  in  the  presence  of  an  accelerating  field  0^ 
behaves  as  in  the  case  of  the  vertical  field  error,  i.e., 
below  threshold  it  exhibits  the  Fresnel  behavior  and 
above  threshold  the  ring  locks  into  the  resonance.  The 
same  is  also  true  during  acceleration  even  in  the  absence 
of  a  vertical  or  toroidal  field  error,  but  in  the  presence  of 
stellarator  fields  of  periodicity  m,  when  the  quasiequili- 
brium  position  of  the  ring  is  off  the  magnetic  axis  of  the 
strong  focusing  system  and  ihe  resonance  mode  is  equal 
to  m.  The  results  of  the  studies  with  axial  field  error  will 
be  reported  in  a  future  publication.  In  the  absence  of  a 
toroidal  magnetic  field,  nonlinear  effects  associated  with 
the  crossing  of  resonances  in  synchrotrons  have  been  con¬ 
sidered  previously  by  DePackh  [6]. 

The  preceding  discussion  is  based  on  the  assumption 
that  the  space  charge  is  low  and  the  strong  focusing  field 
is  zero.  In  addition  to  introducing  new  characteristic 
modes,  the  strong  focusing  field  makes  the  expression  for 
the  regular  cyclotron  mode  more  complicated  [7].  How¬ 
ever,  it  can  be  shown  that  for  the  parameters  of  the  NRL 
device  and  provided  /  » 1,  the  strong  focusing  has  only  a 
minor  effect  on  the  cyclotron  resonance.  This  conclusion 
is  supported  by  extensive  computer  calculations. 

This  paper  is  organized  as  follows.  The  theoretical 
model  is  formulated  in  Sec.  II.  Section  III  A  contains  ex¬ 
amples  of  both  the  Fresnel  and  lock-in  state  from  the  ex¬ 
act  equations  of  motion.  The  slow  equations  of  motion 
are  derived  in  Sec.  Ill  B,  and  the  simplified  slow  equa¬ 
tions  of  motion  with  linearized  detuning  factor  are  given 
in  Sec.  Ill  C.  The  asymptotic  behavior  in  the  Fresnel  and 
lock-in  state  as  well  as  the  appropriate  initial  conditions 
to  be  used  in  the  subsequent  sections  are  presented  in  Sec. 
IV,  while  Sec.  V  contains  a  discussion  of  possible  ways  to 
cross  a  resonance  without  locking  into  it.  Resonance  dia¬ 
grams  for  nonzero  initial  perpendicular  velocity  are 
displayed  in  Sec.  VI,  and  the  multiple  crossing  of  the 
same  resonance  is  demonstrated  in  Sec.  VII.  Finally,  Sec. 
VIII  contains  the  summary  and  conclusions. 

II.  MODEL  WITH  A  VERTICAL-HELD  ERROR 

Imperfections  in  the  coils  that  generate  the  betatron 
field  could  result  in  a  field  error.  A  typical  example  is 
given  in  Fig.  1,  which  shows  the  VF  error  per  kiloampere 
of  the  current  circulating  in  the  coils  as  a  function  of  the 
toroidal  angle,  at  r  =  100  cm.  This  error  is  due  to  a  small 
straight  section  in  each  coil  that  generates  the  betatron 
field  in  the  vicinity  of  the  power  feeds.  Table  I  provides 
the  Fourier  decomposition  of  the  error.  The  values  are 
the  actual  error  amplitudes  (in  G)  for  a  toroidal  field 
Rflf,  —  4650  G  as  each  particular  resonance  is  reached  dur¬ 
ing  acceleration. 

Obviously,  near  a  resonance,  only  the  mode  associated 
with  that  resonance  is  acting  on  the  ring.  The  contribu- 


FIG.  1.  Vertical  field  error  per  kiloampere  of  current  in  the 
coils  that  generate  the  betatron  field,  as  a  function  of  the 
toroidal  angle  6,  at  r  =  100  cm,  z  =0.0  cm. 


tion  of  all  the  other  modes,  being  far  away  from  the  reso¬ 
nance,  averages  out  to  zero  due  to  their  fast  oscillatory 
behavior.  Therefore  the  VF  error  will  be  expressed  in 
terms  of  the  particular  /  mode  associated  with  the  reso¬ 
nance  under  study,  i.e., 


- sin(/0-(-0(j)  , 

'■q 


6R,  —  SBjq 


1  +K^ 


sin(/6-l-0o)  , 


(la) 

(lb) 


TABLE  1.  Fourier  decomposition  of  VF  error.  Actual  values 
of  the  VF  error  at  each  /  mode  for  B^  =  46S0G  are  listed. _ 


Fourier 

mode 

Fourier 

amplitude 

cos  (Fourier 
amplitude) 

sin  (Fourier 
amplitude) 

1 

108.2698 

81.6151 

•’1.1430 

2 

1.2152 

0.3749 

1.1559 

3 

8.8142 

-6.5489 

5.8994 

4 

0.9985 

0.0905 

0.9944 

5 

1.0148 

-0.3835 

-0.9395 

6 

0.5915 

-0.5915 

0.0000 

7 

0.3350 

-0.0856 

0.3238 

8 

0.3188 

0.0884 

-0.3063 

9 

0.3226 

-0.2428 

-0.2124 

10 

0.0759 

0.0640 

0.0407 

11 

0.1301 

0.0868 

-0.0970 

12 

0.0243 

-0.0243 

0.0000 

13 

0.0373 

0.0264 

0.0264 

14 

0.0112 

0.0088 

-0.0069 

15 

0.0062 

-0.0045 

0.0042 

16 

0.0034 

0.0008 

0.0033 

17 

0.0010 

-0.0010 

-0.0002 

18 

0.0012 

-0.0012 

0.0000 

19 

0.0002 

-0.0002 

0.0000 

20 

0.0000 

00000 

0.0000 
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8Bg  =  bB^Q—c(x{l6+6o)  . 
''o 


(Ic) 


In  Eqs.  (1),  6B,q,  9q,  and  are  the  amplitude,  the 
phase,  and  the  gradient  in  the  radial  direction  of  the  / 
mode  and  Tq  is  the  major  radius  around  which  the 
analysis  is  carried  out.  It  should  be  noticed  that  Eqs. 
(la)-(lc)  satisfy  Maxwell’s  equations  to  first  order  in 
toroidal  correlations. 

The  magnetic  fields  acting  on  the  ring  in  the  MBA  are 
the  betatron,  the  toroidal,  and  the  stellarator  field.  For 
simplicity,  the  stellarator  field  is  omitted  in  the  present 
analysis.  At  high  energies  its  contribution  to  the 
confinement  of  the  ring  is  diminished  and  abundant  com- 
puter  runs  have  shown  that  the  main  results  presented 
here  are  not  altered  in  the  presence  of  the  stellarator  field 
except  for  resonances  l  =  km,  where  m  is  the  field  period 
of  the  stellarator  and  fc  =  1,2,  . .  .  .  In  the  analysis,  the 
betatron  field  :s  approximated  by 


B 


(b). 


B 


zO 


(2a) 

(2b) 


where  n  is  the  field  index  and  B^q  the  field  on  the  minor 
axis.  Similarly,  the  jroidal  magnetic  field  is  given  by 


Bg=Bgo^  ,  (3) 

where  Bg^  is  the  field  at  r  =  Tq. 

Since  the  VF  error  is  a  sinusoidal  function  of  the 
toroidal  angle  6,  it  is  convenient  to  express  the  equations 
of  motion  in  terms  of  the  indeoendent  variable  9  rather 
than  time.  In  their  transformed  state,  the  equations  of 
motion  become 


r'+ 


rPe  '  ”  2  !  +  /> 

,ro(n,+/nj 


r 


rP(f 


{l+P)^+\  +  P  ,  (4) 


where  the  complex  variable  ^=[(r  — ro)  +  /z]/ro. 

=  d^/d6,  is  the  complex  conjugate  of  and  fl  is 
the  cyclotron  frequency,  i.e.,  11=  \elB  /me.  The  position 
of  the  ring  centroid  is  given  by 

P^- — ^  =  +  ^  (5ai 

'•o 

e=--=f(^-r)-  (5b) 


The  quantity  ypg  can  be  expressed  as 


rPe^ 


rl-L 


1  +  - 


1 


(1+F)' 


■rr 


\n 


(6) 


where  y  is  the  relativistic  factor  and  its  derivative  with 
respect  to  9  is  equal  to 


rPe 


(y^6)'  = 


(y/5e)' 


>•0^,  1 

rP^  2/ 


(r-r) 


(1+P) 


2  1  +  F 


(7) 


The  accelerating  electric  field  is  given  by 

£g=~~roB,o>  (®^ 

c 

where  B^  is  the  time  derivative  of  the  betatron  field  at 
r  =  ro.  Using  the  rate  of  change  of  time  that  is  given  by 
t’  =  {rQ/cpg){  l+P),  it  can  be  shown  that 


2 

ftzo(l+^) 


and 


(1+P)  . 


(9) 


(10) 


In  the  equations  given  above  fl,  and  fl,  include  both  the 
betatron  field  and  the  VF  error.  Notice  that  Eq.  (4)  is  a 
second-order  nonlinear  differential  equation  of  the  com¬ 
plex  quantity  In  terms  of  the  position  of  the  ring 
centroid  is  given  by  Eqs.  (5),  while  its  normalized  velocity 
components  v,,v^  are 


Pr  +  iP.= 


rPe 

r 


i+p- 


(11) 


Obviously,  the  exact  set  of  nonlinear  equations  given 
above  is  very  complicated  and  difficult  to  handle.  How¬ 
ever,  for  a  ring  that  has  a  bounce  motion  with  a  small 
amplitude,  i.e.,  for  small  mismatch,  and  for  large  /-mode 
values,  it  is  easy  to  show  that  near  the  resonance 
l^i  =  |f'|//  and  therefore  |^i  is  much  smaller  than  if'  . 
In  this  case,  it  is  appropriate  to  linearize  Eq.  (4)  with 
respect  to  f  and  f*.  Under  such  conditions,  Eq.  (4) 
simplifies  to 


(12) 


where 


/,=A  + 


:i-f  A)(f'+f')+/7»f(f  +  f*)- 


(rPe^\ 


y'-|(i-nC)(fr+f*r’)-T<rr+^r') 


-A:,f*-)-fC(f4f)-6C[l+f  f(14-A'K,)i*  + j(r  +  r')r]s'n(/^+^o)  - 


^yPeh 


14f'f' 


11/2 


(I. 3a) 
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(r^e)ic  ’ 

{13b) 

^  0^*0 

irPehc  ’ 

(13c) 

roSft,o 

6C  = -  . 

(rPehc 

(13d) 

2  ’ 


A=l-C  , 

(131) 

(13g) 

/:j  =  (l-n)C-l  . 

(13h) 

When  the  mismatch  term  A  is  small,  then  C  ss  1  and 
K2~\~n.  Since  the  term  in  / ,  will  be  treated  as  a 
perturbation  in  the  derivation  of  the  slow  equations  of 
motion,  must  be  small.  Therefore  we  assume  that  n  is 
close  to  Also,  since  we  consider  only  values  of 
||;'|  <0.2,  before  the  resonance  is  crossed,  terms  proper- 
tional  to  are  omitted  when  compared  to  one.  The 

simplified  Eq.  (12)  is  complemented  bv  the  derivatives 
with  respect  to  0  of  y  and  raSl^Q/c,  which  are  given  by 


Results  from  the  numerical  integration  of  Eq.  (12)  and 
the  exact  Eq.  (4)  are  in  very  good  agreement,  provided 
the  initial  values  of  |^1  <0.30X  10”^  and  |^1<0,2. 
Therefore  Eq.  (12)  will  be  used  in  the  derivation  of  the 
slow  equations  by  averaging  out  the  fast  cyclotron  oscil¬ 
latory  motion. 


Ul.  DYNAMIC  BEHAVIOR  NEAR  A  RESONANCE: 

FRESNEL  AND  LOCK-IN  REGIMES 

A.  Numerical  results  from  the  exact  eqnatioas 

The  results  from  the  numerical  integration  of  the  exact 
equations  of  motion,  i.e.,  Eq.  (4),  as  a  resonance  is  crossed 
during  acceleration  indicate  the  following  consistent  be¬ 
havior:  there  is  a  threshold  value  of  the  VF  error  ampli¬ 
tude  below  which  the  perpendicular  velocity  of 
the  gyrating  electrons  increases  by  a  certain  amount  as 
the  resonance  is  crossed  and  then  it  remains  relatively 
constant,  after  the  resonance  has  been  crossed.  Above 
the  threshold  value,  keeps  increasing  with  time  while 
remains,  on  the  average,  constant,  and  the  resonance 
is  never  crossed,  i.e.,  the  ring  is  locked  into  the  reso¬ 
nance.  A  typical  example  of  this  behavior  is  shown  in 
Fig.  2  for  the  parameters  listed  in  Table  II.  The  thresh¬ 
old  value  of  bB^  is  between  0.19  and  0.195  G.  Figures 
2(a)  and  2(c)  show  Py  versus  time  below  and  above  thresh¬ 
old,  while  Figs.  2(b)  and  2(d)  show  the  corresponding  yPg 
versus  time. 

It  is  rather  difficult  to  find  from  the  exact  equations  of 
motion  the  source  of  the  dynamic  behavior  shown  in  Fig. 
2.  This  behavior  can  be  explained  by  the  slow  equations 
of  motion  derived  in  Sec.  Ill  B. 

To  derive  the  slow  equations  of  motion  the  instantane- 


18.8 - 

CD  ' 

<3Q.  la^ 

18.0- 
0  "  ^  2 


'"i  ^WW  i 

181:  ' 


FIG.  2.  Pi  and  ypg  vs  time  obuined  from 
the  exact  equations  of  motion  [Eq.  (4)],  in  the 
Fresnel  and  lock-in  regimes  and  close  to  the 
threshold,  for  the  parameters  in  Table  II. 
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TABLE  II.  Parnneters  of  the  runs  shown  in  Figs.  2-6. 


Parameter 

Value 

Torus  major  radius  Tq 

l(X)  cm 

Toroidal  magnetic  field  B^o 

2771  G 

Vertical  magnetic  field 

305  G 

Field  index  n 

0.5 

Rate  of  change  of  vertical  field 

2  G/fisec 

ReMnance  mode  / 

9 

Amplitude  of  VF  error 

0.190,0.195  G 

Constant  phase  of  VF  error  6q 

0.0 

Gradient  of  VF  error 

0.0 

Initial  normalized  toroidal  momentum 

17.922 

Initial  normalized  vertical  velocity  Rj 

0.0 

Initial  phase  of  vertical  velocity 

0.0 

Initial  radial  displacement  r  — rg 

0.0  cm 

Initial  vertical  displacement  z 

0.0  cm 

Integration  time  tf 

8  ^lsec 

ous  position  of  the  particle  ^  is  decomposed  into  two 
components.  The  first  is  associated  with  the  bounce 
motion  and  the  second  with  the  cyclotron  or  fast  motion, 
i.e.,  ^=^6+^c.  where  I*  =|4”cxp( 

=  |;^/’exp(  — /Vf  0),  and  v^,  and  are  the  two  charac¬ 

teristic  frequencies  of  the  system.  Specifically,  V4=v_ 


and 


-i(;e+e„) 


(16b) 


It  is  apparent  from  Eqs.  (16)  that  it  is  possible  to  extract 
the  bounce  and  (slow)  cyclotron  motion  in  the  rotating 
frame  by  inverting  Eqs.  (16).  Specifically, 


b 

1 

2 

(15)  v,-Vfc 

Vc 

(17a) 


The  complex  amplitudes  and  are,  in  general, 
slowly  varying  quantities,  provided  that  the  perturbation 
is  not  very  large.  Since  these  two  amplitudes  vary  slowly, 
the  corresponding  derivatives  with  respect  to  0  are  equal 
to  +  provided 

that  =r/.  The  amplitude  is  a  slowly  varying  quanti¬ 
ty  because  has  been  expressed  in  a  frame  that  rotates 
with  angular  velocity  /c/rg.  As  a  consequence  of 
I  I  « 1^'!,  the  position  of  the  cyclotron  mode  is  given 
by  I;,  =  -  ( 1  //V, ) F;‘’exp[  - » ( / 0  -b  0o )  ]. 

From  the  previous  discussion,  the  complex  position  ^ 
and  its  derivative  can  be  expressed  as  follows: 

^  -  ( 1  /«•  V, )  V  "  (1 6a) 


and 


K!*'= 


(tV(,^-t-|’')e 


iae+e^) 


(17b) 


Figures  3(a)  and  3(b)  show  the  real  and  imaginary  com¬ 
ponents  of  00yy'  for  the  parameters  listed  in  Table  11. 
For  these  same  parameters,  the  actual  orbit  in  the  r-z 
plane,  during  the  first  2  fzsec,  is  given  in  Fig.  4(a),  while 
Fig.  4(b)  gives  the  bounce  motion  associated  with  the  ac¬ 
tual  orbit,  as  computed  from  £q.  (17a). 

B.  Slow  equations  of  motion 

The  method  that  is  used  in  this  section  to  obtain  the 
equations  of  motion  of  the  slowly  varying  quantities  asso- 


-0  08  -0  04  0  0.04  0.08  0.12  -0.24  -0.16  -0  08 


FIG.  3.  Normalized  perpendicular  velocity 
in  velocity  space  and  in  the  rotating  frame  |Eq 
(17b)]  obtained  from  the  exact  equations  of 
motion  in  the  Fresnel  and  lock-in  regimes  (for 
the  same  parameters  as  in  Fig.  21. 


0 


0.08 


-40 


2048 


D.  DIALETIS,  S.  J.  MARSH,  AND  C.  A.  KAPETANAKOS 


47 


FIG.  4.  (a)  Exact  orbit  and  (b)  bounce 
motion  orbit  obtained  from  the  exact  equa¬ 
tions  of  motion  during  the  first  2  ^sec,  and  for 
the  same  parameters  as  in  Figs.  2  and  3 


ciated  with  the  cyclotron  and  bounce  motions  has  been 
devised  by  Bogoliubov  and  Mitropolsky  (8)  and  is  a  per¬ 
turbation  approach  to  the  problem.  Before  this  method 
can  be  applied,  the  simplified  Eq.  (12)  should  be  written 

in  a  suitable  form.  For  this  purpose  we  express  in 

terms  of  the  quantities  as  follows: 

+  — J".  .  (18a) 

v+ 

,  (18b) 


where 


i>Q  — 


'’o^eo 


Co  = 


(19a) 

(19b) 

(19c) 


^10~^0  T  * 


1 

(y^fiio 


y2_i 


1/2 


(19d) 

(19c) 


In  Eq.  (19e),  Kq  is  the  zero  order,  slowly  varying  quan¬ 
tity  associated  with  [see  Eq.  (26a)].  In  addition,  are 
definitely  slowly  varying  quantities,  since  they  depend  on 

The  equations  of  motion  of  and  are 


v:+iv^y,- 


v+  -v_ 


(20a) 


/  +  6/  +  ivL^^ 


(20b) 


where 


/  =  Ao+7<l+Ao)<^'  +  r')r-SCoSin(/O+0o)  . 
8/  =  A-Ao  + 


-<(h-ho)+j(A-Ao)(r+r’)+'H($+r)-7-^y 

(yPe'i 


-i(i-nC)(^r+rr'i-T»rr+^r') 


-8C[^  +  (14-A^)r-)-l(^'-)-^*')f']sin(/0  +  ^o>'(SC-6Co)sin(/0  +  0o> 


^0~  1  ^0  ' 


8Co 


^rPe^cf  ' 


r 


and  already  been  defined. 

To  implement  the  perturbation  theory,  we  introduce 
the  parameter  e,  which  indicates  the  relative  smallness  of 
the  various  terms.  In  terms  of  e,  Eqs.  (20)  can  be  written 
as  follows: 

F;-fi/F^=e[/^-((v^-/)FJ  ,  (23a) 


where 


v+  -V.. 


(21a) 


(21b) 

(22a) 

(22b) 


(23b) 


(24a) 
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h=-- 


■[/  +  6/]  . 


(24b) 


Assuming  that  the  field  error  |Aol«l, 

and  since  /  » |  Vq\  and  v+  »/,  the  right-hand  side  of  Eqs. 
(23)  is  at  least  of  order  6.  In  addition,  since  |v_|  «  1,  the 
term  in  Eq.  (23b)  is  of  order  e.  The  terms  that  are 

proportional  to  the  derivatives  of  have  been  omitted 
because  they  are  of  higher  order. 

Furthermore,  for  high  /-mode  numbers  ||;|  «|^|  and 
in  light  of  Eqs.  (18),  in  terms  of  €,  ^  and  ^  can  be  written 
as  follows: 


where  slowly  varying  quantities  associated 

with  respectively.  Then,  the  slowly  varying  quan¬ 

tities  A  1,^2,  ■  ■  ■  and  B^,B2,  ■  ■  ■  are  determined  by  in¬ 
serting  the  equations  above  into  Eqs.  (23)  and  eliminating 
the  terms  that  lead  to  secular  terms  for  each  order  of  the 
perturbation  parameter  e. 

The  computation  is  carried  out  to  second  order  in  e 
and  it  is  tedious  but  straightforward  and  will  not  be  given 
here.  To  that  order,  the  quantities  F,  are  equal  to 


1 


v_  il 


(25a) 

(25b) 


The  perturbation  method  is  applied  on  Eqs.  (23)  and  (25) 
by  setting 

=  •••  .  (26a) 

•  •  •  .  (26b) 


^1  = 


v+— v_  21 


Ao-(-i(l+Ao)lKoP+|6Coe'' 


j(l-t-Ao)^"^ 


(27a) 


(27b) 


Also,  to  second  order  in  e,  the  slow  equations  for  Vq  and 
|:oare 


(26c) 


(26d) 
_ ) 


yo+nv^-nvo=f^^' . 

and 

io+/v_io=/;’' , 

where 


(28a) 


(28b) 


An  + 


''-(^0 


(29ai 


+  V^(bVl-l^'oi')-Ao(l+2^o-Ao<l  +  |lFoP))]Fo 


2 


l+^0+(l+/fK))^0+4l'"ol'-4^'o- 


6Cn 


8v.  l  +  |Kol^  ° 


and 


/r=- 


(/’o'"-  i 


^20 


nCn 


|Kol'  +  (Ao+|Koi')/;o 


^0* 


SCq 

2v4 


^10^0 

2  l  +  |Tol' 


iyo+yv+f^rc^o- 


8Co 


(2%) 


The  quantities  for  and  K20  are  given  in  Eqs.  (J3e) 
and  (13h)  with  replaced  by  iFgP  and  C  replaced  by 
Cq.  The  slow  Eqs.  (28)  and  (29)  should  be  complemented 
by  the  derivatives  of  y  and  rgHjo/c.  Keeping  only  the 
zero  order  nonoscillatory  parts  in  Eqs.  (14a)  and  (14b),  we 
obtain 


(30) 


and 


'■oh; 


20 


n 


lO 


(y^e)o 


1  + 


(31) 


'■0 

2 

^0+^0* 

c 

^20 

+ 

Solutions  of  the  slow  equations  are  shown  in  Fig.  5  for 
the  same  input  parameters  (see  Table  11)  used  in  the  solu¬ 
tion  of  the  exact  Eq.  (4).  Figures  5(a)  and  5(c)  show  the 
zero  order  (y/3e)ol  Fgl /y  versus  time,  below  and 
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FIG.  5.  01  and  y0^  vs  lime  obtained  from 
the  slow  equations  of  motion  [Eqs.  (28)  and 
(29)],  in  the  Fresnel  and  lock-in  regimes  and 
close  to  the  threshold,  for  the  parameters  in 
Table  II. 


Time  (/isec)  Time  (/xsec) 


Fresnel  Regime  Lock— in  Regime 


02468  02468 

Time  (^sec)  Time  (/xsec) 


FIG.  6.  Detuning  factor  vs  time  and  nor¬ 
malized  perpendicular  velocity  in  velocity 
space  obtained  from  the  slow  equations  of 
motion  in  the  Fresnel  and  lock-in  regimes  (for 
the  same  parameters  as  in  Fig.  5). 
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above  threshold,  while  Figs.  5(b)  and  S(d)  show  the  corre¬ 
sponding  zero  order  iY^g)o  versus  time.  Figure  6  shows 
the  detuning  factor  w—v^—l  and  the  real  and  imaginary 
parts  of  below  and  above  threshold.  By 

comparing  Figs.  5  and  6  to  Figs.  2  and  3  it  is  apparent 
that  the  slow  equations  predict  accurately  and  y/3^ 
when  the  initial  perpendicular  velocity  is  zero.  Above 
threshold  the  detuning  factor  locks  into  the  resonance, 
with  only  a  small  variation  around  zero.  A  plot  of  is 
not  given,  since  it  is  exactly  the  same  as  that  in  Fig.  4(b). 
Obviously,  the  numerical  integration  of  the  slow  equa¬ 
tions  is  much  faster  than  that  of  the  exact  £q.  (4)  oi  (12), 
since  the  oscillatory  terms  have  been  averaged  out. 

C.  Slow  equations  of  motion 
with  linearized  detuning  factor 

The  slow  equations  of  Sec.  Ill  B  are  still  too  lengthy  to 
provide  insight  into  the  origin  of  the  threshold  behavior. 
When  Kol//B,o«  1*  the  last  term  in  Eq.  (29b)  can 
be  neglected  and  the  equation  for  the  bounce  motion 
simplifies  to  the  equation 


from  Eqs.  (19e)  and  (34)  that 

(yPa)o=[y^-l-|yol^]'^-  (35) 


A  straightforward  transformation  of  Eq.  (33)  from  the 
Kq  to  the  Uq  variable  leads  to  the  slow  equation  for  (/q, 
namely 


U'^  +  i 


1 


(36) 


where 

„  _  1 


(37) 


Equation  (36)  indicates  that  the  detuning  factor 
w  =  v^—l  should  be  redefined  by  adding  to  it  the  small 
correction  term 


Co 


^0 


1  It/pp 
^  (y^a)o 


/2v.. 


lo+'v_ 


-fj- 


l  +  l^(boV-  +  Ao+l^'oP) 

A  10 


nCo 


^0 


10 


^20 - + 


SO 


v_ 


=  '-I^('^o+t1^'oP)  .  (32) 

A  10 


where  we  have  used  the  fact  that  |v_|  «v+,  and  the  re¬ 
lation  v_v  +  —  ^10'  Equation  (32)  indicates  that  the 

equilibrium  position  of  the  bounce  motion  is  proportional 
to  the  generalized  mismatch  Ao-l-||  FqI*.  Without  excep¬ 
tion  the  computer  runs  have  shown  that  during  the  cross¬ 
ing  of  the  resonance  the  equilibrium  position  of  the 
bounce  motion  hardly  changes  when  the  gradient  factor 
K^=0.  That  is  not  the  case  when  K^¥^0.  Therefore,  if 
both  the  initial  mismatch  and  the  initial  bounce  position 
^0  are  zero,  the  latter  quantity  remains  very  small.  Under 
these  conditions,  the  slow  equation  for  Vq  simplifies  to 
the  equation 


Vk  +  i 


v^-/+^(Ao-j|Fol') 


8C 


/  yy  y _ 2.(  1  -f  '  1 i  y2  \ 

irPe)l  °  2 


(33) 


where  terms  proportional  to  Ag,  |Kol\  and  AqIKoP  or 
higher  have  been  omitted. 

To  simplify  further  Eq.  (33),  it  is  convenient  to  intro¬ 
duce  the  quantity  Uq  by  means  of  the  relation 


Vq^ 


Uo 


[y'-i 


(34) 


Since,  to  zero  order,  Vq  is  equal  to 
[{(3,+iP,)/P0]e%pli(l9  +  6o)],  Uq  is,  to  the  same  order, 
equal  to  yiP^  +  iP,  )exp[i(/0-l-0o)],  and  it  is  easy  to  show 


Above  threshold,  the  small  correction  term  is  important 
because  it  compensates  for  the  small  time-dependent  con¬ 
tribution  of  /Tio  in  v+  [cf.  Eq.  (19a)].  If  the  small  correc¬ 
tion  term  is  omitted  in  either  Eq.  (36)  or  (33),  then,  above 
threshold,  the  average  value  of  (yPe)Q  does  not  remain 
constant  but  increases  with  6.  This  is  contrary  to  the 
solution  of  either  Eqs.  (28)  and  (29)  or  the  exact  (4)  or 
(12). 

When  the  mismatch  Aq + y  I  Kq  I  ^ » 0  and  the  initial  ^q  is 
also  zero,  then  ^q^O  and  the  relativistic  factor  y  is  equal 
to  y  =  yo  +  y'^  [cf.  Eq.  (30)],  where  Xq  y  *be  ini¬ 
tial  values  of  y  and  the  acceleration  rate  [assumed  to  be  a 
constant,  and  given  by  Eq.  (30)  with  respectively. 

Assuming  that  y'9«yo  and  HUqI^—  IUqqI^I  «},  where 
Uqq  is  the  initial  value  of  Uq,  the  detuning  factor  can  be 
linearized  with  respect  to  y'0and  |Uool^  namely 


iu  =  iuo-a0-l- j6[|t/oP-|Uool^]  , 

where 


“^o  =  '^+o~/  + 


'00 


2v+o 


Aflo 


1  ic/, 


00' 


^  (y^eW 


a  =  v 


YoY 

""'(y^elro  ’ 


6  =  v.^ 


I 


v..  =  -2 


(yPe)M 

hoo 


l-f 


^00 

2 


^00 

2 


-i-a: 


100 


TTT 


(38) 


(39a) 


(39b) 


(39c) 


(39d) 


^yPb^oo  [Yo  (  i^ool  1  • 


(39e) 


and  v+o.hoo*Coo' Aoo-^ioo  the  initial  values  of 

v+,6o,Co,  Ao,Af  10,  respectively.  The  dependence  of  K  ^Q 
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FIG.  7.  0^  and  ypg  vs  time  obtained  from  the  slow  equations 
of  motion  without  acceleration,  for  the  parameters  in  Table  III. 


on  6  and  on  iy^e^o  neglected,  since  v.^  is  weakly 

dependent  on  K  ig  when  bo^l»l.  Also,  since 
Q-l.  and  b~l/(y0g)lo,  the 
small  correction  term  that  is  added  to  the  detuning  factor 
in  Eq.  (36)  provides  such  a  small  contribution  to  6,  of  or¬ 
der  l//^  that  it  has  been  neglected.  Therefore  the  slow 
equation  for  Uq  with  linearized  detuning  factor  is 

C/i  +  i(u;o-a0+iS(|l/ol^-|C/ool')]t/o=-<Bo  • 

(40) 

This  equation  predicts  the  same  temporal  behavior  for 
and  ypg  as  the  original  slow  Eqs.  (28)  and  (29).  The 
cause  of  this  behavior  is  the  nonlinear  term  proportional 
to  shown  shortly  that  when 

8=0,  the  solution  is  the  Fresnel  integral,  which  is  bound¬ 
ed).  We  conclude  that  the  existence  of  the  Fresnel  and 
lock-in  regimes  in  the  exact  Eq.  (4)  or  in  Eq.  (12)  is  due  to 
the  dependence  of  the  term  r^ngf^/yff^  upon  Y0g  and  the 
fact  that  y0g  depends  nonlinearly  on  the  perpendicular 
velocity. 


It  is  interesting  to  examine  the  case  when  there  is  no 
acceleration  (a=0),  the  particles  are  at  exact  resonance 
(ii>o~0)  and  the  initial  perpendicular  velocity  is  zero 
(I  t/ool^=0).  Then  Eq.  (40)  becomes 


l/|,+/j8|t/oPC/o=-/5o  . 


(41) 


and  it  can  be  solved  exactly.  For  that  purpose,  define 
A  =  |l/oP.  Then,  separating  the  real  and  imaginary  parts 
of  Eq.  (41),  it  is  easy  to  show  that 


Vor=- 


A^  , 
A’  , 


(42a) 

(42b) 


and  from  the  definition  of  A  it  follows  that 

I  1/2 


A'  = 


—  —  A  ^ 

4JIo  A 


(43) 


The  exact  solution  of  this  differential  equation  is 

/)J  A  dx  _ 

0  [xd-xUx^+x  +  D]'^  ^  ’ 


(44) 


where  >.  =  (6/8Bo)''^^.  The  integral  can  also  be  expressed 
in  terms  of  the  elliptic  integral  of  the  first  kind  [9] 
F{(p,k),  i.e.. 


F 

2  arccot  — ~ 

\~k^A 

1/2  - ^ 

v'2-v^3 

31/4 

k^A 

’  2 

(45) 


For  very  small  or  very  large  values  of  the  argument  of  the 
arccot,  the  approximate  solution  is 

ii/J 


A  =(ro^p^= 


85. 


ian  fi6 


v^  +  tan^pS 


(46) 


where  p  =  3''^^A.5o-  The  solution  of  Eq,  (43)  is  bounded. 
When  p5  «1,  A^{Bq6)^  and  is  proportional  to  6, 
which  is  the  usual  linear  secular  solution.  But  the  non¬ 
linear  dependence  of  yPg  on  forces  the  solution  to  be 
bounded.  Figure  7  shows  and  yPe  versus  time  for  the 
parameters  listed  in  Table  III,  by  integrating  either  the 


TABLE  111.  Parameters  of  the  run  shown  in  Fig.  7. _ 

Parameter  Value 


Torus  major  radius 

100  cm 

Toroidal  magnetic  field  Bgo 

2771  G 

Initial  generalized  mismatch  Ao-f0.5|Fol^ 

0.0 

Field  index  n 

0.5 

Rate  of  change  of  vertical  field 

0.0  G/psec 

Resonance  mode  / 

9 

Amplitude  of  VF  error  6B,o 

0.3  G 

Constant  phase  of  VF  error  do 

0.0 

Gradient  of  VF  error 

0.0 

Initial  normalized  toroidal  momentum  yPg 

18.1753 

Initial  normalized  vertical  velocity  5i 

0.0 

Initial  phase  of  vertical  velocity  ^0 

0.0 

Initial  normalized  radial  bounce  displacement 

0.0 

Initial  normalized  vertical  bounce  displacement  ^oi 

0.0 

Integration  time  tf 

8  pscc 

47 


DYNAMIC  BEHAVIOR  OF  AN  ELECTRON  RING  CLOSE  TO  A  .  .  . 


2053 


slow  £q.  (40)  or  (36),  or  £qs.  (28)  and  (29).  Integration  of 
the  exact  Eq.  (4)  for  the  same  parameters  also  gives  iden¬ 
tical  results.  The  peak  value  of  /9|  and  the  period 
(ir/^)fQ/c  as  computed  from  the  approximate  Eq.  (46) 
arc  0.0757  and  1.247  /rsec,  respectively,  in  good  agree¬ 
ment  with  Fig.  7. 

In  conclusion,  our  analysis  indicates  that  there  is  a 
completely  different  behavior  near  a  resonance  with  or 
without  acceleration.  Even  at  exact  resonance,  0^  is  al¬ 
ways  bounded  without  acceleration,  but  it  is  unbounded 
with  acceleration  and  in  the  presence  of  a  large  field  er¬ 
ror. 


rv.  ASYMPTOTIC  BEHAVIOR 
IN  THE  FRESNEL  AND  LOCK-IN  REGIMES 


In  order  to  study  the  asymptotic  behavior  of  Eq.  (40),  it 
is  convenient  to  transform  it  into  dimensionless  form. 
For  this  purpose,  when  a  >  0,  we  introduce  the  quantities 
a  =  v^d,  ao=WQV^,  and 

Oo  =  {b/2a''^)Uo. 

Then,  Eq.  (40)  becomes 

^+i[ao-<r-^\Oo\^-\Doo?]Oo=-ii,  (47) 

where  Uq^  is  not  equal  to  the  initial  value  of  Dq  and  will 
be  defined  shortly.  Let  be  the  threshold  value  of  f. 
When  i.e.,  when  .Bo8/(2a^''^)  the  change 

in  Oq  associated  with  the  crossing  of  the  resonance  is 
small,  i.e.,  — l^ool^l  «  1,  and  Eq.  (47)  simplifies  to 


dOp 

da 


i{ap  —  a)0p=  —ii . 


(48) 


initial  condition  of  Eq.  (SO),  and  also  of  Eq.  (40),  by 
means  of  the  initial  variable  Ugp,  i.e.. 


Uop  =  0, 


00® 


liWg/ltt) 


Wp 


(51) 


where  Dpp  =  (6/2a''^)~'0pp  and  Uqq  is  the  value  in  Eq. 
(47)  around  which  the  detuning  factor  is  linearized.  Then 
Eq.  (50)  becomes 


.<»  /2a) 


00 


J-  ^0  f 


(52) 


and  is  independent  of  Tq,  but  it  does  depend  on  the  initial 
phase  of  Equation  (52)  indicates  that  as  6  tends  to 
minus  infinity,  1 1  I  ^ool^l  tends  to  zero.  Therefore 
IC/oqI  is  the  appropriate  constant  around  which  the  de¬ 
tuning  factor  should  be  linearized,  and  Wp,a,h  become 
functions  of  I/qo  (and  not  of  Because  |  t/ol  tends  to 
l&ool  ^  ~  0°  we  shall  call  £7qq  the  asymptotic  initial 

value  of  Vq.  When  0»1  (or  |u)l/Va»l),  the  asymp¬ 
totic  value  of  I  I/q  I 

Wol-lUoo-iBpv'vTad+i)]  .  (53) 


Therefore,  when  f/oo^O,  the  final  value  of  IUqI  could 
be  smaller  than  its  asymptotic  initial  value  ll/ooi  When 
Uq^,  —  0,  the  asymptotic  value  of  0^  is 


/3i(e»l)  =  v^2tT/a 


Bp 

n+r'^  ’ 


(54) 


where  \y'6\  «Yp.  Since  the  width  of  the  Fresnel  integral 
is  Aa  =  2fl''^^,  we  conclude  that  the  time  At  it  takes  for 
the  resonance  to  be  crossed  is 


The  crossing  of  the  resonance  occurs  at  a=ap,  i.e.,  at 
^p—Wq/u.  The  solution  of  this  equation  is 


(55) 


Hi/2)icr-a^)^ 


0, 


00® 


-1(1/2 1<7^ 


(49) 


i.e.,  the  particular  solution  can  be  written  in  terms  of  the 
Fresnel  integral.  When  cro»l,  the  asymptotic  expres¬ 
sion  of  the  Fresnel  integral  leads  to  the  approximate  solu¬ 
tion 


B„ 


Wn 


~HujQ/2a) 


B 


1 


.  "0  r  - 

Va  J- 


(50) 


where  the  solution  has  been  transformed  to  the  original 
variables  and  iv^Wg—aff.  Uq  is  very  sensitive  to  the 
phase  u>o/2a,  when  l/pQ+Bg/tVp^^O.  Since  a«I,  a 
small  change  in  the  initial  time  interval  Tq  =  (  Tq  /c  )wq  /a 
from  the  resonance  [e.g.,  a  small  change  in  the  initial 
value  of  (y^glol  will  cause  a  very  different  behavior  of 
I  t/fl  I  versus  time.  To  alleviate  this  problem,  we  define  the 


Therefore  the  final  0^  and  At  are  inversely  proportional 
to  the  square  root  of  the  acceleration  rate  in  the  Fresnel 
regime,  provided  i  «  ^thr-  Figure  8  shows  0^^  and  {yBe^p 
versus  time  for  the  parameters  of  Table  II,  except  that 
6Bjq=0.  1  G,  i.e.,  far  away  from  the  threshold.  Notice 
the  smooth  variation  of  during  the  first  microsecond 
and  compare  to  the  oscillatory  behavior  that  occurs  in 
Fig.  5(a).  The  difference  is  due  to  the  choice  of  the  initial 
condition.  In  Fig.  8(a),  while  in  Fig.  5(a), 

Ugg  =  0.  When  w>0  and  \w\/'^a»l,  i.e.,  well  before 
the  resonance  crossing,  the  asymptotic  value  of  Ug,  as 
given  by  Eq.  (52),  is 


I!  ~-r! 
t/n  '-'00® 


(56) 


When  Opp—0,  we  have  jf/ol  =  — Bq/IuiI,  which  explains 
the  smooth  variation  of  0°^  in  Fig.  8(a>. 

We  have  shown  above  that  if  the  initial  condition, 
given  by  Eq.  (51),  is  chosen  the  solution  of  Eq.  (40),  in  ei¬ 
ther  the  Fresnel  or  the  lock-in  regime,  is  independent  of 
the  injtial  time  interval  Tg  from  the  resonance,  when 
Wg/'^a»l.  This  is  no  longer  true  if  the  initial  condi¬ 
tion  given  by  Eq.  (51)  is  chosen  for  the  solution  of  Eq. 
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FIG.  8.  Pi  and  yPe  vs  tin's  from  the  slow  equations  of  motion 
in  the  Fresnel  regime  far  away  from  the  threshold,  i.e.,  for  the 
parameters  of  Table  II,  except  that  5B,o“0-  •  G  [ths  initial  con 
dition  is  given  by  Eq.  (51),  with  i/oo=0]. 

(36),  which  has  a  nonlinear  detuning  factor  w.  The  non¬ 
linear  dependence  of  to  on  0  causes  the  solution  to  depend 
on  Tq.  To  second  order  in  6,  the  detuning  factor  is  given 
by 

u;  =  u)o-a0  +  iaj0‘-i-i8[|L^o!^-|t^ool^l  •  <57) 

where  lOg.a.S  have  already  been  defined,  C?oo 
redefined  shortly,  and 

_  ,  (y^e)oo  nr' 

a-,=  3 - ; —  - —a  .  (58) 

rl 

When  w  /^a»\,  i.e.,  well  before  the  resonance  cross¬ 
ing,  the  term  proportional  to  5  can  be  neglected  in  Eq. 
(57),  and  the  asymptotic  solution  of  Eq.  (40)  is 


where 


Inverting  Eq.  (59b),  we  get  (since  026^ « 1 ) 


ii'  -Wq  a.  , 

0= - 1 - r  (  W  -  Wq  ?■  , 

a  la^ 


[  r  ^  1 

—  /  /  w  d  6 

1  •'o 

a9+\a20^  ■ 

U)o 

and,  therefore. 


/•  d 

w  d6=  I  4- 

n 


CjICq  Wq 

6a^  2a 


then  Eq.  (S9a)  becomes 

w 

Notice  that  Uf,  depends  on  ujq  and,  therefore,  on  Tq. 
However,  if  a2WQ/a^ «  I,  then  the  dependence  is  weak 
and  the  solution  of  Eq.  (36)  with  the  initial  condition 
given  by  Eq.  (62)  is  also  weakly  dependent  on  Tq.  Equa¬ 
tion  (63)  indicates  that  as  0  tends  to  minus  infinity, 
llf^ol^“l^ool^l  tends  to  zero,  and,  therefore,  Ugg  is  the 
appropriate  parameter  around  which  the  detuning  factor 
should  be  expanded.  Also,  iUo,a,a2,5  are  functions  of 
Oqq.  As  before,  we  call  Oqo  the  asymptotic  initial  value 
of  t/o-  From  this  point  on,  the  solutions  of  Eqs.  (28)  and 
(29)  are  obtained  using  the  initial  condition  computed 
from  Eq.  (62).  If  such  an  approach  is  not  followed,  each 
small  change  in  the  initial  parameters,  e.g.,  iyPe^oo’-Beo’ 
etc.,  would  produce  a  different  Ug  which  could  be  either 
in  the  Fresnel  or  the  lock-in  regime,  to  the  extent  that 
one  might  get  the  impression  that  the  behavior  is  ran¬ 
dom. 

An  analytic  expression  of  on  0,  above  threshold, 
can  be  easily  obtained  assuming  that  iyPe^o  remains  al¬ 
most  constant.  This  assumption  is  justified  by  the  numer¬ 
ical  solution  of  either  Eqs.  (28)  and  (29)  or  Eq.  ',36).  Re¬ 
sults  are  shown  in  Fig.  9(b)  for  the  same  parameters  listed 
in  Table  II,  except  that  8Bjo=0.4  G.  In  contrast  to  the 
solutions  of  Eqs.  (28)  and  (29)  or  Eq.  (36),  Eq.  (40)  pre¬ 
dicts  that  iy0g)o  initially  behaves  in  the  same  manner,  but 
gradually  increases  as  9  (or  time)  increases  to  very  large 
values.  The  exact  Eq.  (4)  or  (12)  gives  also  an  average 
yPg  that  varies  as  in  Fig.  9(b).  Since  ypg  remains  con¬ 
stant,  on  the  average,  after  the  resonance  has  been 
reached,  we  conclude  that 

(To+y’^)^“  1  “  i<^o!^~consi  .  (64) 

If  the  initial  value  t/oo=0,  then  at  0  =  9g  =  Wg/a,  we  have 
I  t/o  I  ==  0  and  T  =  To  T  ^O’  Fq.  (64)  leads,  to  zero  or- 


0.2 


Cti 

o 

18.18 

18.10 

H--y(lO0-4ti'') 

a 

18  02 

2a 

(61) 

17.94 

Instead  of  Eq.  (51),  where  a2~0>  'f  we  redefine  the  initial 
condition  by  the  relation 


^00 


i(  1  ^  a^wQ/ba^  )uQ/2a  "0 


Time  (/xsec) 


FIG.  9.  Asymptotic  Ijehavjor  of  0^  and  yP^  above  the  thresh¬ 
old  (lock-in  regime)  for  (he  parameters  of  Table  II,  except  that 
6fi,o=0.4G 


4” 
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der  in  to  the  expression 


1 4-  — (0+0o) 


1+^0 

Yo 


Yo 


(65) 


When  y'0«Yo>  fii  is  proportional  to  the  square  root  of 
time,  which  is  demonstrated  in  Fig.  9(a)  and  also  is 
proportional  to  the  square  root  of  the  acceleration  rate 
Y'-  In  contrast,  below  threshold,  it  has  been  shown  that 
Py  is  inversely  proportional  to  the  square  root  of  y'. 


V.  DYNAMIC  STABii^iZATION 
AND  THRESHOLD  LAW 

There  are  at  least  three  possible  ways  to  cross  a  cyclo¬ 
tron  resonance  without  locking  into  it:  (a)  by  reducing  or 
eliminating  the  held  errors,  (b)  by  accelerating  fast 
through  the  resonance,  and  (c)  by  adding  a  small  time- 
dependent  toroidal  field  that  provides  dynamic  stabiliza¬ 
tion.  In  effect,  it  will  be  shown  that  the  dynamic  stabili¬ 
zation  is  equivalent  to  increasing  the  acceleration  rate. 
The  effectiveness  of  the  second  stabilizing  mechanism  is 
discussed  at  the  end  of  this  section.  Here,  we  analyze  the 
effectiveness  of  the  stabilizing  toroidal  field,  which  is  as¬ 
sumed  to  be  sinusoidal  with  amplitude  63^^,  and  period  r 
which  is  much  larger  than  the  time  At  is  takes  to  cross 
the  resonance.  The  total  toroidal  field  is  described  by 

i?9o(t)=£^+6JS9oSin - —  ,  (66) 

T 


FIG.  10.  Threshold  law  in  the  presence  of  a  stabilizing  time- 
dependent  toroidal  magnetic  field,  i.e.,  6B«=&B«osin(2n’i/r) 
and  zero  initial  perpendicular  velocity  (&(io=0). 


rate,  the  larger  becomes  the  threshold  value.  Therefore  a 
small  time-dependent  toroidal  field  with  negative  time 
derivative  does  provide  dynamic  stabilization. 

The  threshold  law  is  obtained  by  determining  numeri¬ 
cally  the  threshold  value  of  i  in  the  dimensionless  Eq. 
(47)  in  the  special  case  when  Uoo—0.  Then  the  threshold 
law  is  expressed  by  the  relation 


where  is  the  time  delay.  Inserting  tq.  (66)  into  Eqs. 
(13b}  and  (19b)  [with  the  initial  value  £^(0)  in  bool  and 
linearizing  the  detuning  factor  in  Eq.  (36),  the  parameters 
0,02  become 


o  =  V 


-n 


T 


(67a) 


,  (y/5e)w 

Intj 

02  = 

Yo 

(y/3flW 

T 

.  (69) 

Figure  10  gives  the  threshold  values  of  as  a  function 
of  6£jq  for  different  acceleration  rates,  when  £35  =  3771 
G,  t=50  /xsec,  and  /  =9.  The  bei  fleial  effect  of  a  large 
acceleration  rate  or  a  small  time-dependent  toroidal  field 
becomes  obvious  from  this  figure.  Each  line  in  the  figure 
separates  the  Fresnel  from  the  lock-in  regimes  for  each 
different  value  of  the  acceleration  rate. 


where 


(67b)  VI.  DYNAMIC  BEHAVIOR 

WITH  INITIAL  PERPENDICULAR  VELOCITY’ 


Toy' 

_  In  rp  Yo 

r  c  ’ 


(68a) 

(68b) 


and  u)o,6,v+|  are  given  by  Eqs.  (39a),  (39c),  and  (39d). 
The  expansion  of  the  detuning  factor  in  Eq.  (36)  is  done 
around  and  not  around  the  initial  value  C/qo.  so  that 
all  the  parameters  given  above  depend  on  Ugo-  For  sim¬ 
plicity,  let  us  assume  that  tj=0.  If  the  toroidal  field  de¬ 
creases  as  the  resonance  is  being  crossed,  i.e.,  if  85 <0, 
then  according  to  Eq.  (67a)  a  becomes  larger.  This  is 
equivalent  to  replacing  y’  in  a  with  a  larger  effective 
value.  It  is  shown  shortly  that  the  larger  the  acceleration 


When  the  asymptotic  initial  value  of  Uq  is  not  zero,  the 
display  of  the  dynamic  behavior  close  to  the  resonance 
becomes  more  difficult,  because  the  dynamics  depend  on 
the  initial  phase  as  well  as  on  the  initial  amplitude  of  Dq. 
The  results  are  conveniently  presented  as  contour  plots 
(i.e.,  the  resonance  diagrams)  of  the  final  P.  =  \Uq\  /y  in 
the  0^\O),{pQ  plane,  where  0y°HO)  and  <Pq  are  the  ampli¬ 
tude  and  phase  of  the  asymptotic  initial  value  ^co^Yo- 
The  rest  of  the  parameters,  including  the  interval  of  in¬ 
tegration  over  6  or  time,  are  kept  constant.  Since  the 
final  py  becomes  large  in  the  lock-in  regime  while  it  is 
bounded  in  the  Fresnel  regime,  the  contours  are  very 
dense  at  the  boundary  between  the  two  regions. 

Figure  1 1  shows  the  two  regimes  for  the  parameters 
listed  in  Table  IV.  Although  the  results  shown  in  Fig.  1 1 
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Table  IV.  Parameters  of  the  nm  shown  in  Figs.  1 1  and  12. 


Parameter 

Value 

Torus  major  radius  Tq 

100  cm 

Toroidal  magnetic  field  5^ 

2771  G 

Initial  generalized  mismatch  Ao-t-O.SlFgP 

0.0 

Field  index  n 

0.5 

Rate  of  change  of  vertical  field 

2,4  G//iscc 

Resonance  mode  / 

9 

Amplitude  of  VF  error  65^ 

0.195,0.3  G 

Constant  phase  of  VF  error  5o 

0.0 

Gradient  of  VF  error 

0.0 

Amplitude  of  stabilizing  toroidal  field  65«o 

0,-200  G 

Period  of  stabilizing  toroidal  field  r 

70  psec 

Time  delay  of  stabilizing  toroidal  field 

0.0  fliCC 

Initial  normalized  toroidal  momentum 

17.922 

Initial  normalized  radial  bounce  displacement 

0.0,0.003,-0.003 

Initial  normalized  vertical  bounce  displacement 

0.0 

Integration  time  t( 

8  /isec 

have  been  obtained  by  integrating  Eqs.  (28)  and  (29),  Eq. 
(36)  or  (40)  gives  the  same  results,  provided  that  there  is 
negligible  bounce  motion.  In  Fig.  11(a),  has  been 
chosen  equal  to  0.195  G,  which  is  very  close  to  the 
threshold  value  when  ^/’*(0).  By  increasing  the  field  er¬ 
ror  amplitude  from  0.195  to  0.3  G,  the  lock-in  regime 
dominates  for  small  initial  0°\O)  for  the  entire  range  of 


initial  phase  angles  as  shown  in  Fig.  11(b).  By  turning  on 
the  stabilizing  field,  the  Fresnel  region  increases  at  the  ex¬ 
pense  of  the  lock-in  region.  Results  arc  shown  in  Fig. 
1 1(c)  for  dB^=0. 195  G  and  f>BgQ—  — 2(X)  G.  This  figure 
should  be  compared  with  Fig.  11(a)  that  has  the  same 
field  error  amplitude  but  not  stabilizing  field.  Finally,  in 
Fig.  11(d)  the  acceleration  rate  has  been  increased  from 


FIG,  11.  Resonance  diagrams,  i.e.,  contour  plots  of  the  final  Bi.  when  the  initial  amplitude  and  phase  of  the  normalized  perpendic¬ 
ular  velocity  Dgn/yo  are  not  zero,  for  the  parameters  in  Table  IV,  and  in  (al  6B,o=0-195  G,  8Bgo=Q  O  G,  B;p  =  l  G/jusec;  (b) 
6B,o  =  0  3  G,  6B«,=0  0.  B^  =  1  G/^xsec;  (c)  6fl,o=0195  G,  8fl8o=-200  G,  t  =  70  fxsec,  5,0=2  G//isec;  (d)  65, i,  =  0.195  G. 
65^0  =0.0  G,  5,0  =  4  G//xscc. 
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=  2  to  4  G//isec,  while  =  195  G  and 

Again,  comparison  of  Fig.  11(d)  with  Fig.  11(a)  shows 
that  the  Fresnel  region  has  increased  at  the  expense  of 
the  lock-in  region.  When  the  electrons  in  the  ring  are 
uniformly  distributed  over  the  initial  phase  angle,  the  res¬ 
onance  diagrams  give,  for  each  initial  the  percen¬ 

tage  of  the  ring  that  crc.,.-;es  the  resonance  and  the  per¬ 
centage  that  locks  into  i,. 

In  all  the  results  presented' so  far,  the  equilibrium  posi¬ 
tion  of  the  bounce  motion  was  located  at  the  origin  of  the 
coordinate  system  [i.e.,  Aqo+ I t/ool^/2(y^g)^— 0]  and 
the  amplitude  of  the  bounce  motion  was  selected 
negligibly  small  (since  the  initial  *he  equi¬ 

librium  position  is  chosen  off  the  origin  [i.e., 
I and  the  amplitude  of  the 
bounce  motion  is  negligibly  small,  by  judicially  choosing 
|;oo,  the  resonance  diagrams  remain  the  same  for  the  same 
set  of  parameters.  When  Aoo+ |C/ool^/2(y/3g)^9^0,  and 
in  the  special  case  of  a  field  index  «  =  the  parameters 
y'  and  rg/c  in  Eqs.  (68a)  and  (68b)  should  be  replaced  by 
(l+A,o)y'  and  (l  +  Aiglrg/c,  respectively,  where 
!  ^ool^/2(y^5)^]//L  log.  Also,  the  initial 
bounce  position  must  be  equal  to  ^oo“‘^io*  to  have 
bounce  motion  with  negligible  amplitude.  Finally,  the 
term 


0  100  200  300 


Vq  (deg) 

FIG.  12  Resonance  diagrams  for  the  same  parameters  as  in 
Fig.  1 1(a),  except  that  the  radius  of  the  bounce  motion  is  0.3  cm 
in  both  la)  and  in  (b)  and  in  (a!  the  initial  position  of  the  bounce 
motion,  I.e.,  ^^*,  =  (0.3  cm,  0.0  cm);  in  (b)  £no~(~0.3  cm,  0.0 
cm). 


1  1  / 

(v_  26g?o)  |Ao+  2  j/^iO 

should  be  added  to  the  detuning  factor  in  Eq.  (36),  and, 
therefore,  the  term  (v_o“26go^oo)A|o  should  be  added  to 
Wq  in  Eq.  (39a),  where  v_o  is  the  initial  value  of  the 
bounce  frequency. 

When  there  is  a  small  bounce  motion  superimposed  to 
the  cyclotron  motion  there  is  a  modest  change  of  the  res¬ 
onance  diagrams.  This  becomes  apparent  by  comparing 
Fig.  11(a)  with  Fig.  12.  These  figures  have  been  obtained 
from  Eqs.  (28)  and  (29)  with  the  same  parameters,  except 
for  the  initial  bounce  position.  In  Fig.  1 1(a),  the  initial 
position  of  the  bounce  motion  is  at  the  origin,  while  in 
Fig.  12(a),  it  is  at  (0.30  cm,  0.0  cm)  and  in  Fig.  12(b)  it  is 
at  (  —  0.30  cm,  0.0  cm).  In  the  latter  two  cases,  the  ampli¬ 
tude  of  the  bounce  motion  is  0.30  cm. 

Vn.  MULTIPLE  CROSSING  OF  A  RESONANCE 

When  a  small  toroidal  field  with  sinusoidal  time  depen¬ 
dence  is  added  to  the  main  toroidal  field  the  detun¬ 
ing  factor  may  become  zero  more  than  once,  as  time 
evolves,  for  the  same  resonance  mode  /.  This  is  defined  as 
a  multiple  crossing  of  the  resonance  mode  /.  Results  are 
shown  in  Fig.  1 3.  These  results  have  been  obtained  from 
Eqs.  (28)  and  (29)  for  the  parameters  listed  in  Table  V  and 
S5^o=0.35  G.  Figure  13(a)  shows  that  the  detuning  fac¬ 
tor  w  becomes  zero  five  times  and  a  Fresnel  jump  occurs 
each  time  the  resonance  mode  /  is  crossed.  It  is  not 
necessary  that  all  the  crossings  have  a  Fresnel  jump.  A 
lock  into  the  resonance  could  occur  at  some  crossing  if 
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FIG.  13.  V,  — yBfi,  and  0^  vs  time  dunng  a  multiple  cross¬ 
ing  of  the  resonance  /  =  8,  obtained  from  the  slow  equations  of 
motion  [Eqs.  (28)  and  (29)j,  for  the  parameters  of  Table  V,  and 
68,0=0.35  G. 
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TABLE  V.  Parameters  of  the  run  shown  in  Fig.  7. 


Parameter 

Value 

Torus  major  radius 

100  cm 

Toroidal  magnetic  field 

2771  G 

Initial  generalized  mismatch  Ao 

0.0 

Field  index  n 

0.5 

Rate  of  change  of  vertical  field  B,o 

0.8  G/fiscc 

Resonance  mode  / 

8 

Amplitude  of  VF  error 

0.35,1.0  G 

(Constant  phase  of  VF  error  8o 

0.0 

Gradient  of  VF  error  K,q 

0.0 

Amplitude  of  stabilizing  toroidal  field  SBgo 

-200  G 

Period  of  stabilizing  toroidal  field  r 

20  fiscc 

Time  delay  of  stabilizing  toroidal  field 

0.0  fiscc 

Initial  normalized  toroidal  momentum  yPb 

19.875 

Initial  normalized  vertical  velocity  0^ 

0.0 

Initial  phase  of  vertical  velocity  q>o 

0.0 

Initial  normalized  radial  bounce  displacement  ^q. 

0.0 

Initial  normalized  vertical  bounce  displacement 

0.0 

Integration  time 

80  fiscc 

the  asymptotic  initial  velocity  is  in  the  lock-in  regime  of 
the  resonance  diagram.  As  a  rule,  the  detuning  factor 
follows  the  time-dependent  toroidal  field.  Figure  13(b) 
shows  that  when  all  the  crossings  are  in  the  Fresnel  re¬ 
gime,  Y^e  follows  Y  between  crossings,  while  Fig.  13(c) 
shows  that  there  are  five  Fresnel  jumps  in  0^.  In  contrast, 
when  the  ring  locks  into  the  resonance,  then  YPe  follows 
the  time-dependent  toroidal  field  while  w  remains  very 
small.  This  is  shown  in  Fig.  14,  obtained  by  integrating 
Eqs.  (28)  and  (29),  for  the  parameters  listed  in  Table  V 
and  6Bjo“  10  G.  According  to  Fig.  14(a)  the  first  cross¬ 


ing  occurs  at  1  fiscc  and  the  ring  locks  into  the  resonance 
up  to  approximately  10  fiscc.  As  long  as  it  remains 
locked,  u)  wO,  while  yPb  follows  the  time  varying  toroidal 
field  [Fig.  14(b)].  Figure  14(c)  shows  that  just  before  10 
/usee,  decreases  to  zero,  while  Y^e  increases  at  a  faster 
rate  than  y.  Since  at  this  time 
=  1 -(y^e)^«0,  a  continuously  rising  y^b 

would  require  a  negative  (y/Si)^,  which  is  an  unphysical 
situation.  Thus  the  ring  unlocks  from  the  resonance. 
The  same  cycle  is  repeated  up  to  45  fiscc,  when  /?|,  due  to 
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FIG.  14  —  I,  and  0  vs  time  during  a  multiple  cross¬ 

ing  of  the  resonance  /  =  8  under  the  same  conditions  as  in  Fig, 
13,  except  that  LOG 
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FIG.  15.  V,  —  /,  yPg,  and  0.  vs  time  during  a  multiple  cross¬ 
ing  of  the  same  resonance  /  =  8,  under  the  same  conditions  as  in 
Fig.  13,  but  obtained  by  integrating  the  exact  equations  of 
motion  [Eq.  (4)]. 
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the  rise  of  its  average  value,  cannot  become  zero.  After 
this  time,  yPe  remains  locked  to  the  sinusoidal  variation 
of  the  field,  while  w  remains  almost  zero.  Therefore  the 
next  resonance  mode  /  + 1  is  never  reached,  in  spite  of 
the  fact  that  y  keeps  increasing.  The  energy  is 
transferred  to  yP^^  rather  than  yp^.  The  exact  Eq.  (4)  or 
(S)  or  the  nonlinear  slow  £q.  (36),  give  identical  results  to 
those  obtained  from  Eqs.  (28)  and  (29)  which  arc  shown 
in  Fig.  14.  However,  this  is  not  the  case  with  Fig.  13. 
Results  from  the  integration  of  the  exact  Eq.  (4)  for  the 
parameters  listed  in  Table  V  and  6R,o=0.35  G  arc  shown 
in  Fig.  15.  The  detuning  factor  and  ypg  in  Fig.  15(a)  and 
15(b)  are  very  similar  to  the  results  shown  in  Figs.  13(a) 
and  13(b).  However,  the  transverse  velocity  in  Fig.  15(c) 
is  similar  to  that  of  Fig.  13(c)  until  the  third  crossing 
occurs  at  19  /zsec.  The  difference  is  due  to  the  fact  that 
the  slow  Eqs.  (28)  and  (29)  are  approximate,  and  the 
phase  of  P^  computed  from  the  exact  and  the  slow  equa¬ 
tions  becomes  gradually  different  for  long  times.  There¬ 
fore,  for  long  times,  when  many  resonances  are  to  be 
crossed  or  multiple  crossing  occurs,  the  exact  Eq.  (4) 
should  be  used.  The  slow  equations  provide  a  valuable 
insight  in  the  dynamic  behavior  close  to  a  particular  reso¬ 
nance,  but  are  not  reliable  over  long  periods  of  time. 

VIII.  SUMMARY  AND  CONCLUSIONS 

Magnetic-field  errors  excite  resonances  that  the  elec¬ 
tron  ring  must  cross  during  acceleration.  In  the  presence 
of  a  VF  error,  there  is  a  threshold  value  of  the  field  error 
amplitude  that  separates  two  distinct  regimes.  Below 
threshold  (Fresnel  regime)  and  for  zero  initial  t/oo,  the 
perpendicular  velocity  increases  by  a  finite  amount  as  the 
resonance  is  crossed.  The  increase  as  well  as  the  time  it 
takes  to  cross  the  resonance  are  inversely  proportional  to 
the  square  root  of  the  acceleration  rate.  Above  threshold 
(lock-in  regime),  the  perpendicular  velocity  is  proportion¬ 
al  to  the  square  root  of  the  acceleration  rate  and  increases 
with  the  square  root  of  time,  while  yPg  remains  on  the 
average  constant,  and  the  detuning  factor  remains  ex¬ 
tremely  small.  Therefore  the  ring  locks  into  the  reso¬ 
nance.  The  dynamic  behavior  without  acceleration  is  en¬ 
tirely  different  even  at  exact  resonance.  The  perpendicu¬ 
lar  velocity  is  proportional  to  time  initially,  but,  due  to 
the  nonlinearities  in  the  equations  of  motion,  it  reaches  a 
maximum,  then  it  decreases  to  zero  and  repeats  periodi¬ 
cally  the  same  cycle.  Therefore  it  is  bounded. 

The  threshold  is  predicted  by  the  slow  equations  of 
motion  that  have  been  derived  by  averaging  out  the  fast 
cyclotron  motion.  The  origin  is  the  nonlinear  depen¬ 
dence  of  yPg  on  the  perpendicular  velocity  and  the  fact 
that  the  cyclotron  frequency  is  inversely  proportional  to 
yPg.  By  the  appropriate  choice  of  the  initial  conditions  it 
was  shown  that  the  solutions  of  the  slow  or  the  exact 
equations  of  motion  could  be  made  weakly  dependent  on 
the  initial  time  interval  from  the  resonance.  Possible 
ways  to  increase  the  threshold  have  been  discussed.  It 
has  been  shown,  that  the  threshold  value  of  the  VF  error 
amplitude  is  proportional  to  the  ^  power  of  the  accelera¬ 
tion  rate.  Dynamic  stabilization,  t.e.,  the  addition  of  a 
small  time-dependent  field  to  the  main  toroidal  field,  also 


provides  an  effective  increment  to  the  acceleration  rate,  if 
it  has  a  negative  time  derivative,  and,  therefore,  increases 
the  threshold.  When  the  initial  perpendicular  velocity  is 
not  zero,  the  dynamic  behavior  has  been  presented  by 
means  of  the  resonance  diagrams.  These  diagrams  pre¬ 
dict  that  a  small  bounce  motion  has  only  a  modest  effect 
on  the  Fresnel  and  lock-in  regimes.  Finally,  the  multiple 
crossing  of  the  same  resonance  has  been  analyzed  in  the 
presence  of  dynamic  stabilization  and  we  have  concluded 
that  for  long  periods  of  time  the  exact  equations  of 
motion  should  be  used. 

Following  the  successful  demonstration  of  acceleration 
in  the  NRL  device,  a  concerted  effort  was  made  to  locate 
and  eliminate  or  reduce  the  various  field  disturbances 
that  may  excite  the  cyclotron  resonance.  Reduction  in 
many  of  these  field  errors,  together  with  the  operation  of 
higher  toroidal  and  strong  focusing  fields  led  to  beam  en¬ 
ergies  in  excess  of  20  MeV,  while  the  trapped  current  was 
above  1  kA. 

In  addition,  three  different  cyclotron  resonance  stabili¬ 
zation  techniques  were  tested  in  the  NRL  modified  batat- 
ron:  enhancement  of  the  acceleration  rate,  dynamic  sta¬ 
bilization  or  tune  jumping,  and  avoidance  of  the  reso¬ 
nance. 

The  damage  done  to  the  beam  at  each  resonance  de¬ 
pends  on  the  speed  with  which  the  resonance  is  crossed. 
By  enhancing  the  acceleration  rate  the  resonance  is 
crossed  faster  and  thus  the  damage  inflicted  to  the  beam 
is  reduced.  To  achieve  higher  acceleration  rate,  the  verti¬ 
cal  field  coils  were  divided  into  two  halves  with  midplane 
symmetry  and  powered  in  parallel.  The  experimental  re¬ 
sults  show  a  striking  reduction  of  the  beam  losses  at 
/=12,  11,  and  10,  when  the  acceleration  rate  increased 
from  0.69  to  1.93  G//isec. 

The  crossing  of  the  resonance  can  be  also  speeded  up 
by  modulating  the  toroidal  magnetic  field  with  a  rapidly 
varying  ripple.  This  is  the  dynamic  stabilization  or  tunc 
jumping  technique  and  requires  a  carefully  tailored  pulse 
to  be  effective  over  many  resonances.  These  results  have 
been  reported  [10]  previously  and  in  general  they  are  in 
agreement  with  the  predictions  of  Sec.  V  and  extensive 
computer  calculations. 

It  is  apparent  from  the  resonance  condition  that  when 

cyclotron  resonance  is  not 
excited.  To  test  this  prediction,  a  linearly  rising  torinda! 
field  ramp  ABp  was  superimposed  on  the  main  loroidal 
field.  During  the  rise  time  of  the  ramp  i  -  100  ^isec.  the 
ratio  + AB^/fi,(,~^onstv^  integer.  The  experimental 
results  indicate  that  during  this  lime  penod  the  beam 
losses  are  completely  suppressed.  Although  very  power¬ 
ful,  the  resonance  avoidance  technique  by  keeping  the  ra¬ 
tio  integer  is  not  practical  because  to 

be  effective  over  the  entire  spectrum  of  I  requires  very 
high  toroidal  field.  Among  the  three  stabilization  tech¬ 
niques  tested,  acceleration  of  the  beam  at  a  higher  ac¬ 
celeration  rate  appears  to  have  the  highest  practical  po¬ 
tential. 
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